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Theory of Angular Correlations 
by F. Coester and J.M.Jauch. 
Department of Physics, State University of Iowa, Iowa City, Iowa. 
(30. X. 1952.) 


Abstract. The most general angular correlation function may be written in the 
form W = Tr (€@), where g is the density matrix describing the ensemble of final 
nuclei and reaction products and the “efficiency matrix” € is completely deter- 
mined by the counter arrangement and the nature of the particles to be counted. 
Expansion of € and 9g into sums of irreducible tensors greatly simplifies the evalu- 
ation of the trace. Explicit general rules are derived for the construction of any 
angular correlation function from a few simple building stones. 


I. Introduction. 


Since the work of Myrrs?) and Hamiiron?) on angular distribu- 
tions and angular correlations many papers*® have been written 
which deal with the general theory of angular correlation of nuclear 
radiations. Their main concern is with the derivation of explicit 
formulae for various correlation functions. The probabilities for 
transitions between definite pure states specified by a complete set 
of quantum numbers are calculated, usually by perturbation me- 
thods, and the angular dependence is derived from the transforma- 
tion properties of the relevant matrix elements. In actual experi- 
ments one observes not individual systems in a pure state but a 
statistical ensemble and not all the observables specifying a pure 
state of the system are actually observed. In the final formula one 
therefore sums or averages over unobserved quantum numbers. 

Considerable simplification in the mathematical formalism can be 
achieved by using the density matrix o*) to describe the ensemble 
of systems each consisting of all the products of a nuclear reaction 
or cascade decay®). With 9, the statistical average of any observable 
and the probability for any observable effect can be calculated. We 
are particularly interested in a “‘yes-no experiment” in which a 
counter arrangement detects all systems in a state specified by the 
set of quantum number @ with an efficiency €g; 0 < €g < 1. The 
probability for the answer “‘yes”’ is 


W= Zieq(@lei®. (1) 
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We may consider €, to be the eigenvalues of an “‘efficiency matrix” 
€ and write (1) in the form 


W = Tr(Eo) (2) 


where the trace in (2) may be evaluated in any representation. Most 
convenient is a representation specified by the quantum numbers 
jm of the total angular momentum and other quantum numbers col- 
lectively denoted by A, A’... Bhe evaluation of the trace is facilit- 
ated if we introduce the irreducible tensors of and €; by the de- 
finition 

(jA lon |’A") = (-1)"-™jj’m—m'| kn) (jmA | o|j’m'A’) (8) 


mm 


and 
(GA | ef] 7A’) = DI (—1)-™"(f7’'m—m'| kn) (jmA]€|j’m' A’) (4) 


where the first factor on the right hand side is a Clebsch-Gordan 
coefficient for the addition of angular momenta®). According to their 
definition €; and o; are matrices with respect to 7 and A but not 
with respect to m; they shall be called efficiency tensors and statisti- 
cal tensors respectively. Their knowledge is equivalent to the know- 
ledge of the density matrix and the efficiency matrix, for with eq. 
(A2) of Appendix A there follows from (3) and (4) 


(jmA | g| 7’ m'A’) = D1 (—1—™ (jj’m—m'| kn) (GA | of|7’A’) (5) 


kn 


and 
(jmA |e | j’m'A’) = FF (—1)'-™ (jj’m—m'| kn) GAJER| FA). (6) 


An 


From (5) and (6) follows with (A1) and (A3) 
W=Tr(Eo)= 2) DS’ DAA) Al |7'A)*. (7) 


jj) AA’ kn 
We also use the fact that the density matrix is hermitian. 


The €; will depend explicitly on the geometry of the counter 
arrangement while @; will depend on the geometric characteristics 
of the initial states. If we consider the decay of isotropically distrib- 
uted nuclei the density matrix must be isotropic, that is diagonal in 
m and independent of m. Hence we find by (3) (j4 | of | 74’) = 0 
for k +0. If we consider a nuclear reaction between non polarized 
particles we have cylindrical symmetry about the beam. The den- 
sity matrix is diagonal in m if the direction of the beam is chosen as 
axis of quantization. There follows from (8) (7A | oe; | 7A’) = 0 for 
n +0. 
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(7) is the most general angular correlation function. In order to 
get formulae useful for the interpretation of specific experiments 
we must derive explicit expressions for the oe; and the ¢€?. This 
will be done for the 0? in Sec. IT and for the €7 in See. ITI. 


II. The Density Matrix. 


The density matrix @ in (1) describes the end products of a nuclear 
reaction or a cascade decay. The initial nuclei are described by the 
density matrix (jmA y|@|4’m’Ay’) or by the statistical tensors 
(j Ag | 0%| 7’Ao’). The probability for a transition from a state Q to 
a state P during the time ¢ is for large t equal to t/2z|Rpg|? where 
S = R + 1is the so called Hetsenperc S-Matrix’). If interactions 
with external magnetic fields of macroscopic or atomic origin are 
negligible§) (jmA|S|j’m’ Ao) is invariant under rotations of the coor- 
dinate system. It is therefore diagonal in 7 and m and independent 
of m and may be written in the form (A|S;|A9). The density matrices 
o and @ are related by the equation®) 


Te 


(jmA |o|j’m' A’) = const. SY (jmAg|o|i’m'A))R;(A,Ap)R}(A’, Ad). (8) 
Ay Ao’ 


Consequently we have with (3) 


(jA |ot|7’A’) = const. 3” (jg | 08 |j’A¢) Rs (A, 4) RP (A',,AQ). (9) 
° Ao Ad’ 

In the case of the decay of isotropically distributed metastable 
nuclei (jAy|o7\j' Ao) is different from zero only for one value of j, 
7’ = 7, and one value of Ay, Ap’ = Ap. Ag can therefore be dropped 
from the notation and we have 


ree 1 
) m | Q | ym ) = 27+1 Om m’ (10) 


and hence ae 
(j | of 19) =V2j+1 4,, Sn - (11) 


If the initial nuclei are partially aligned (j|o%\7) will be different 
from zero for 0 <k < 27, —k <n <k. These tensor elements 
describe the state of polarization. 

We proceed to determine 0% for nuclear reactions. 7 is now the 
total angular momentum of the target nucleus plus incoming par- 
ticle which may be a photon or a heavy particle (proton, «-particle, 
neutron, etc.). If a compound nucleus is formed j is the spin of the 
compound nucleus. The spin quantum numbers of the target nucleus 
are j, m,. L is the total angular momentum of the incoming particle 
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with respect to the target nucleus; the quantum numbers are L, M. 
The center of mass of the entire system is supposed to be at rest. 
For particles of definite energy one more quantum number besides 
L and M is needed to specify a state. This will be the parity p*) 
for photons and the orbital angular momentum / for heavy particles. 
We use & as a short notation for both L, p and L, l. Instead of the 
angular momentum quantum numbers the direction of the momen- 
tum 2 = (0, ¢) and the component of the spin in that direction o 
may be used to specify the state of one particle. The transformation 
matrix between the two representations is given by 


(Qo|&M) = pa (Oo | M’) Digy (9, 8, 0)* (12) 


where Dixy (y, 8, x) are an irreducible representation of the rotation 
through the EuErR angles ¢, 8, 1°). For heavy particles with spin s 
(0o|L 1M) is given by 


(13) 


which vanishes unless M =o. For magnetic multipole photons, 
p = (— 1)", one finds 


(00| LM(—1)**4) =—(L100|LM) ‘ES -o|/* oy, tt) 


i= 
(14) vanishes for o = 0 as required by the absence of longitudinal 
photons. For electric multipole photons, p = (— 1)¥, 


(0o|LIi(—1)") = (11 00| LM) eo® C1 (15) 


where the coefficients c,; are to be determined by the requirement 
that the right hand side of (15) vanishes foro =Oand J’ c? =1. 
This gives l=L+1 


(00|LM (—1)) = 6,4 yee (16) 


A beam of particles moving in the direction 2 is described by a 
density matrix 


(wEM |9(Q)| o'&'M’) =) d(o—o')0, -(0)( Qo |EM)*(Qo"| é'M’) (17) 
where ps 


Tr 0 = [00S ecel (w) =1 (18) 


*)p= agn 12+ for magnetic multipoles, p=(— 1)” for electric multipoles. 
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0,0’ describes the polarization of the beam") and the energy distri- 
bution in the beam. To simplify the notation it will be convenient 
to drop the 6-function with respect to the energies and have it 
understood that the matrix elements are functions of the energy. 
Since all matrices discussed in the following will be diagonal in the 
energies this procedure will be used throughout. If the beam is 
practically monoenergetic it is appropriate to treat @ as a fixed 
parameter and normalize the density matrix by 


a Quo = 1. (19) 
According to (3), we have with (A8), (C5) and (C4) 
(&| of (2) | &’) ar p> Goa (@) (00| Fa) (00'| go’) x 
x (— 1)! (LE o—o' |kn’) DE. (, 8,9) . (20) 


If the particles are not polarized o,,, (w) = e(w) 6,, and 


(§ lot | &) = o(@) C, (&, €”) Dig (9, 9, 0) 
where 
(—1)2+8-s Se ee eC ee : 
C.= > aes V(2l+1) (20° +1) (2L+1) (2L’+1) x 
x (1U'00|k0) WLI LL’, ks) (22) 


for heavy particles and 


C, = 7 V@L +1 @L'+1) > x 


x (1+ pp’(—1)*) (—1)¥-! (LL/1—1|k0) (28) 


for photons. The definition of the Racah coefficient W which appears 
in equ. (22) is given in Appendix B. 

The target is described by a set of statistical tensors (j;| o%¢ | j1) 
which vanish for k, + 0 if the nuclei are not aligned. The statistical 
tensors of the entire system can now be constructed as follows 


Ege onl T Si) =D Delete | de) (Flore | &) x 
ke me ko Mo 


 (jj” [Lieje! LL) | Keke) (Keo meg | ken) (24) 


where J’ is the transformation coefficient defined in (B6) of Appen- 
dix B. 
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For heavy particles it is possible to describe the spin and the 
orbital motion of the incoming particles by separate density ma- 
trices and combine the angular momenta according to the scheme: 
j,. +s =I1,1+1U=j. Lis the so called channel spin. 

(11708 | op {FLU Ges) = DD? D Gelort| ae) (sloes| 8) x 


kent ks Ns ko No 


x (ox Q) |’) SY (ID |Pjeie8 8K) | ky kes) x 


KN 
x (kpkegnyns| KN) (jj |P'(LI 1 k) | Kg) (KlkyNng|kn) (25) 


where 


(Ion (2) |) = <2" VAT) AIF) (11 00| ko) Dig, 9,0). (26) 


This representation is particularly useful only if both target and 
beam are unpolarized, that is if 


(Je | ogt | 41) a Dr, 0 Ongo y2 i:+1 
and 
(8 073 | 8) = bj, 5n,0 2 (@) - 


With (B6), (B4) and (C3) one obtains then from (25) and (26) 


(7 11 js | of (2) | 7’ IU 98) = 
vai (25+1) (27’+1) (2141) (2U’+1) (-1)" ,,,, . 
= On yerpe (23,+1)(21+1) 42 (1U'00|k0) x 
40 
2k+1. 








x Wilk, 17’) ya). (29) 
It is, of course, always possible to transform a j L1s representation 
to a 7Ils representation by (B5). 


Ill. The Efficiency Matrix. 


In order to evaluate explicit expressions for the efficiency tensors 
we must specify the quantum numbers denoted by A, A’ in detail. 
Let € stand for the quantum numbers J, p or L, 1 of any emitted 
particle. Different particles in the same procegs are distinguished by 
subscripts 1,2,3.... We denote now the total angular momentum of 
the system by j, and define j, as the angular momentum of the 
nucleus after emission of the particle 1, j, =j, — L,. Similarly jf, = 
Jo — Le. Jas Jo Je... ave the nuclear spins of successive levels in a 
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cascade decay. We specify the quantum numbers A in successive 
steps: ja, A = Ja §1, J», B where B stands for the remaining quan- 
tum numbers, if any. j,, B = jy, &,j,, C. This procedure may be 
repeated until all emitted particles are described explicitly. 

Since the efficiency matrix for the entire system is equal to the 
product of the efficiency matrices of the individual fragments we 
have 


(7A | pa lig = 2) dy (io Ble 95 BY) (Fx | ps (2x) | &) 


kins ky ny 


* (ja da| Loy Ly Ly Ka) | Ko hey) (Hin bey My My | Kea Ma) (30) 


(E| ER | &’) = 2 Di Soo ( @) (00| &a) (00'| &’ 0’) 
x (—1)¥-" (LL'o—o' |kn’) DE,.(g, 8, 0). (31) 


The derivation of (31) is quite analogous to that of (20). @ is the 
particle energy in the laboratory system. If the counter efficiency 
does not depend on the polarization of the emitted particle we have 
Ce’ (a) Ty. €(@) Oa0’ and 


(5 \€p| &’) = € (@) Cy (&, 8") Dig (gp, 8; 0) « (32) 


According to our derivation 2, in equ. (30) is the direction of propag- 
ation of particle 1 in the center of mass reference frame. For high 
energy nuclear reactions this reference frame is not identical with 
the laboratory frame. Q, is then a function of the corresponding 
solid angle ®, in the laboratory system depending on the energies 
of the incoming and outgoing particles w, and @,. If no other par- 
ticles are observed we have 


(Jo B lene 9 j)B . Or, 0 Ono V2 Jo “4 Oi, jp’ Opp : (33) 


On the other hand if more particles are emitted and observed we 
proceed in analogy to (30) 


(jo B\ ERP | 15 BY) = 2 OY eC | Eee 176 ©") (Fa | Ets (22) | £2) * 


kez ke ne 


(jo95|LUe Fe Le Ls kp) | Ke kg) (Ke hig ee | ky Np) « (34) 


If particle 2 is the last particle observed €;2 is given by (38) with 
proper change of subscripts, otherwise we may reduce it further 
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according to the scheme employed in (80) and (34). Q, is the direc- 
tion of the relative velocity of the particle 2 and the nucleus c. This 
is in general not identical with the direction ®, observed in the 
laboratory system. Q, is a function of ®, depending on the energies 
Wo, @, and w». Q, sad 2, are the observed angles only if the recoil 
of the nucleus is negligible at all stages of the process. 


IV. Summary of Results. 


The development of the preceding sections allows us to construct 
explicit angular correlation functions from (7). We use two sets of 
subscripts to distinguish between the various particles and nuclei. 
Subscript a, b,c, .... label respectively the entire system and the 
residual nuclei after emission of particles labeled by subscripts 1, 2, 
3... The target nucleus in a nuclear reaction is labeled by the sub- 
script t and the incident particle by 0. With (7) and (9) we write the 
general correlation function in the form 


W =const. fro SY Pg dD, (ja A| 2 | 5g 4’) x 


Jaiqg AA’ Ay Ao kana 


(jaAg| Of | ja Aq)* Ri, (A, Ao) Ry, (A’s A’). (35) 


/ dw indicates integration over all energy variables. In the case 
of cascade decay of metastable nuclei j, = j,’ and 0,” is given by 
(11) if the nuclei are not aligned; otherwise Ore describes the align- 
ment. In the case of a nuclear reaction 9% is given by (24) with 
(20) or (21). Numerical values for the J”s can be found from (B6) 
and BrepENHARN’s Tables for the Racau coefficients!?). 

Explicit expressions for €¢ are obtained in a series of steps des- 
cribed by eq. (80) ff. The energy variables w; appear in the following 
places: (1) the energy distribution in the beam @,,/(mp9). (2) the 
energy dependence of the counter efficiencies €,,- (w,). (3) the energy 
dependence of the angles 2; introduced through the transforma- 
tion to laboratory angles ®,. (4) the energy dependence of the 
matrix elements R. For sharp levels the energy dependence of the 
matrix elements can be effectively replaced by 6-functions; for 
broad levels some other plausible assumptions must be made. Selec- 
tion rules and plausible assumptions about the matrix elements 
greatly simplify the formulae to be compared with experiment in 
concrete cases. The discussion of these aspects is, however, outside 
the scope of the present paper?). 
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Two examples should suffice to illustrate the form in which corre- 
lation functions appear. Let us first consider the directional corre- 
lation between two successive radiations from a metastable nucleus. 
Assuming that the intermediate level has a definite spin 7, we obtain 
with the rules. described above 


W (Q,, 22) ) = CON he det hee. eg, (1, €;) C,, (Gas 5) x 


kony hi m ke Me 
x Dk (1 9,0) Di 4(@z 2.0) (jada |L(v jo Ly Ly 0) | ky ky) (Hep ky my m, | 00) x 
X (JoJo Mie je Ly Ly ky)| Ok) (O Fea 0 ng | ky Mp) x 
Bj, (1 §2) By, (&, §5)- (36) 


With A(3), C(3) and C(4) there follows?4) 


W(Q,, 2;) = const. a a 2 Os (&,&;) Cy (&, 5) Px(cos 9) x 


&,6)" &6 


oa: re bias 
-. (ja da\L (ojo Ly L, 0) | kk) x 


V2k+1 
X (JoJo | Le je Ly Ly k)|Ok)R i, ($1 §2) By a G1» §2) 


where @ is the angle between the directions 2, and 2. 


As a, second example we take a nuclear reaction with two emerg- 
ing particles. The beam is nonpolarized and the counters are insen- 
sitive to polarization?). 


W(2)2,2,) = const. [0 [6 Oy, [der 0 (ap) € (@,) € (@e) x 


X 6(@y—.@,— @, + AE) OY dy oy Oe, De, de de dy (& &) x 


&,6," 362" €3 63’ jaja’ join ko Os 


x Cy, (&, €;) Cy, ( (§ &5) Die (Po Fy 0) Di (p, 9, 0) Di (P2920) x 
X (jada |L Ged Lo 145 ko) | Oe) ada | Lode La Ly ho) | ka hy) x 
X (kegky Ng My | Kg Mo) (jndp|LGede Le Laks) | Oh2) RF (S051 S240) x 
« Bi (E51 F240) - (38) 
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V. Internal Conversion. 


According to the theory of internal conversion!) the matrix 
element for the ejection of a conversion electron can be written in 
the form 

>» (Qe|a|LMpi)(LMpb| A |a) (39) 


LMp 


where i denotes the initial state of the electron. a and b are nuclear 
states and (LM pb|H|a) is the matrix element for the emission of a 
photon in the transition a > b. « depends only on an electronic 
wave function and universal functions but not on nuclear wave 
functions. Since 


SY (Qo|a|LMpi)* Qo’ |a| L'M’ pi) = 


= 3, 5 (0o|2|LN pi)*(00"|a|L/N’ p's) Diy( 990) Diy (P0)* (40) 


a 


we may apply all previous results to correlations with conversion 
electrons provided we replace (31) by 


(é|€n| é’) a: oe YD fox (00 |a| LN pi) (0o|a|L'N’p'i) x 


x (—1)"—"" (LL' N—N’ |kn’) DE, (p90) . (41) 
If the polarization of the conversion electron is not observed €,,, = 
6.o’- For K-conversion of pure multipole radiation the coefficients 


in (41) have been calculated by Rosr, BrepENHARN and ARFKEN?’). 
The result can be written in the form 


(& | | &) = by (€) Cx (E16) Dio (p 90) (42) 
where C;, is given by (23) and b, is defined by 


425 FE |(0o|a| LN pi)|?(LLN-N|k0)(-1)"-* 
k= _ SS 





(2 L+1) (-1)"*? (LL 1-1|k0) (48) 


The coefficient b, denoted b, by Rosr, BrepENHARN and ARFKEN 
has been evaluated numerically. 3’ in (43) is taken over K-electrons 
only. ; 
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CLeBscu GorDAN Coefficients?8)?9), 


The Ciesscu Gorpan coeffcients for the addition of angular mo- 
menta a+6 =c are denoted by (ab «| cy). «, B, y are the 3-com- 
ponents of a, b,c. (ab «f | cy) vanishes unless «+ B =y and ¢ is 
one of the numbers a+b, a+b—1, ..., |a—b|. The CresscH 
GorpAN coefficients are real and satisfy the orthogonality relations 


& (aba p| cy) (abaB\c'y’) =6,.4,, (A1) 
ap 


D> (aba B\cy) aba’ B’ | ey) = 6,455 78 ,445° (A2) 


They also satisfy the following relations 


and 


{aba B|cy) =(—1)*"°~* (ba Ba| ey) =(—1)**"""(ab—a—B|c—y) (AB) 


and 


(abaB| ey) = |/2¢+* (1) (aca—y|b—) = 


oy (— 1)°** (cb—y B | a—a) ’ 


Appendix B. 


Racau Coefficients?®)!%) and Transformations between Angular 
Momentum Coupling Schemas. 


The Racan coefficient W(abced; ef) is related to the CLEBscH 
GorDAN coefficients by 


W(abed; ef) 6,,.6,, 3 Sea (aba B|ee) x 


apd ep 
x (eded|cy) (bdBd|f¢y) (afag|c’y’). 
It satisfies the symmetry relations 


W(abcd; ef) = W(bade; ef) = W(cdab; ef) = 
= W(acbd; fe) = (—1)’**-*! W(aefd; be) 


and the sum rule 
D> (29+1) (—1)st?tet4 e149 W(abde; eg) W(acdb; fg) = 
g 
= W(abcd; ef). (B3) 
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If any one variable vanishes the value of the Racaun coefficient can 
be written down immediately 


b+e—f 
W(abed; Of) = {— Dav Xa 
ite am VaiuGen ; = 


Let a,b,c be three angular momentum vectors;d =a+6+ ce; 
e=a+b, f=a--c. The transformation matrix between the coup- 
ling schemes abecd and acfbd is given by 


(e|I"(abed)|f) = 2) 2 (ba Bu|ee) (ecey|dd) x 


aBy ep 


x (acay| fy) (bf By |dd) = V(Qe+l) (2f+1) Wabed; ef). (BS) 


Let a, a’ 6b’ be four angular momentum vectors a+ b=c; 
a’+b’=c';c—c’=d; a—a'=e; b—D' =f. The transformation 
matrix between the couling schemes abe a’'b’c’ d and aa’e bb’fd is 
given by 

(cc’ |I'(aa’ bb’ d) pp pp (abaB|cy) (a’'b’a’ B’\c'y 

vy’ &@ 


x (—1)*-* (aa’a—« ee) ( 1)”-* (bb’ B—’ lf) x 
x (—1)"~” (ce’ y—y’ | dd) (efe p| dd) = 





= V(2e+1) (2f+1) (2c+1) (2c +1) J’ (2x+1) We'cef; dx) = 


x W(acbf'; bx) W(b’c' ae; a’ x) . (B6) 


Appendix C. 


The Representations of the Rotation Group. 


We list here some often used formulae involving the representa- 
tions Diy (y, 8, x) of the rotation group. If F4/(Q) is any irredu- 
cible tensor depending on the direction 2 = (#, m) we have the 


relation 
FU(2 = 2) Fi (0) Dinw (y, 8,0)*. (C1) 


This holds in particular for spherical harmonics and with 
Y(0) =6 att (C2 
there follows 1 (0) = Oo ae (C2) 


Diol 02)* =/5 573 YE(0,9) E(x) =P, (cos 8) (C8) 
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where P; (cos #) is a Legendre Polynomial. D4,,,, satisfies the rela- 
tions 


Di (yp Ox)* = (—1) 4 DE ye (9 Ox) (C4) 


and 
Di uw (P9 x) Di u(Y 4) = 
-: = Py (L, L,M,M,| LM) (L,L, M;,M,| LM’) Diy (phx). (C5) 


L MM’ 
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Anisotropie der Dielektrizitaétskonstante des Eises 
von F. Humbel, F. Jona*) und P. Seherrer. 
(27. XI. 1952.) 


Summary. From large single crystals of ice, grown by a method published ear- 
lier®), square plates were cut with the major surfaces either parallel or perpendicular 
to the optical axis. 

The dielectric constants and the loss angles were measured as a function of 
the temperature and of the frequency. The dielectric constant ¢, measured in the 
direction of the c-axis was found to be slightly larger than the dielectric constant 
€1, measured in a direction perpendicular to the c-axis. The relative difference 
(Eo— €4¢)/€ at 1 ke/sec is about 15% at —5° C, about 12% at — 20° C and about 
8% at —40° C. 

The dependency of the dielectric constants ¢, and ¢,, and of the loss angles 6, 
and 6,, on frequency is similar to that found in earlier observations on poly- 
crystalline ice. 


I, Einleitung. 


Das Eis, ein im gesamten Naturgeschehen tiberaus wichtiger 
Stoff, zeigt ein eigenartiges dielektrisches Verhalten, dessen Kennt- 
nis von fundamentalem Interesse ist. Die Dielektrizitaétskonstanten 
als Funktion der Frequenz und der Temperatur zeigen eine aus- 
geprigte Dispersion. Die Kurven gleichen sehr stark denjenigen 
polarer Fliissigkeiten, auch hinsichtlich des Temperatureinflusses. 
Die Frequenz, bei welcher die Dispersion auftritt, nimmt mit sin- 
kender Temperatur ab. Die ersten Untersuchungen veréffentlichte 
ErrerA!) im Jahre 1924, seine Resultate wurden nachtraglich 
durch analoge Messungen bestitigt, die von GRANIER?), WintTScH’), 
Smytx und Hircucock’) und Murpny®) ausgefiihrt wurden. Alle 
diese Untersuchungen beziehen sich auf das dielektrische Verhalten 
einer polykristallinen Eismasse, wobei stillschweigend immer an- 
genommen wurde, dass sich die Eisproben aus einer Vielzahl sta- 
tistisch orientierter Kristallite zusammensetzten. Diese Annahme 
ist nicht immer als zutreffend anzusehen. Die Diskrepanzen 
zwischen den Messergebnissen verschiedener Autoren diirften 
zum Teil damit erklart werden. 


*) Jetzt am Pennsylvania State College, State College, Pennsylvania. 
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Eine andere Fehlerquelle ist in der Anwesenheit von Verunreini- 
gungen verschiedener Art zu suchen. So bemerkte GraniER?), dass 
fiir Eis, hergestellt aus Wasser mit kleinen Verunreinigungen von 
H,S80,, die Dielektrizitatskonstante (DK) bei hohen Temperaturen 
und tiefen Frequenzen sehr stark ansteigt. Spiater stellten auch 
Smytu und Hrrcucock?) fest, dass Eis aus einer 0,0002 molaren 
KCl-Lésung (1 KCl-Ionenpaar auf 278000 Molekiile H,0O, d. h. 
1 g KCl auf 67100 g H,O) eine viel gréssere DK als das reine Eis 
aufweist, und zwar betragt die Zunahme bei — 1° C und 300 Hz 
etwa einen Faktor 2,4. 

Die Untersuchungen von OptatTKA®) beweisen, dass Eis aus nicht 
entgastem Wasser grosse Raumladungen enthalt, waihrend reines, 
entgastes Wasser zu keinen Raumladungen fihrt. Die Raum- 
ladungen kénnen unter Umstinden die statische DK um das 
30fache vergréssern. 

GranrIER?), OpLATKA®) und auch andere Autoren behaupten, dass 
ganz reines Eis eine statische DK von 80—82 wie das Wasser haben 
sollte, obwohl meistens héhere Werte gemessen wurden. Eine kiirz- 
lich erschienene Arbeit von Auty und Coin’), deren Untersuchun- 
gen fast gleichzeitig mit den unsrigen ausgefiihrt wurden, befassi 
sich besonders mit dem Einfluss von Verunreinigungen und von 
inneren Spalten auf die DK des Eises und gibt fiir dessen statische 
DK bei 0° C den Wert 91,5 an. 

Das Eis kristallisiert hexagonal. Es ist deshalb eine Anisotropie 
der DK zu erwarten, indem der in Richtung der c-Achse (optische 
Achse) auftretende Wert nicht notwendigerweise gleich sein muss 
dem Werte der DK in einer beliebigen Richtung senkrecht zur c- 
Achse. Die vorliegende Untersuchung setzt sich zum Ziel, die Grés- 
senordnung dieser Anisotropie zu bestimmen. Zu diesem Zwecke 
mussten verschieden orientierte Eiseinkristalle hergestellt und auf 
die Frequenz- und Temperaturabhangigkeit ihrer DK untersucht 
werden. 


II. Apparatur und Messverfahren. 


Die Messung von Kapazitaét und Phasenwinkel des Eiskonden- 
sators erfolgte mit Hilfe einer Briickenschaltung nach Fig. 1. 

C, ist ein geeichter Luftdrehkondensator, verinderlich zwischen 
26 und 109 pF. Die Kapazitit dieses Kondensators ist auf + 0,05pF 
ablesbar. C, ist ein Luftkondensator von 100,0 pF, zu welchem 
jeweils verschiedene, geeichte Zusatzkapazititen C%, parallel ge- 
schaltet werden kénnen. Dadurch erhalt man die Erweiterung des 
Messbereiches fiir C,. C, und Cy sind zwei gleiche Kapazititen von 
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je 1000 pF, bzw. nach Umschalten von je 0,1 wF. Auf diese Weise 
erhilt man die Méglichkeit der Erweiterung des Messbereiches von 
R,. Der Widerstand R, ist ein Dekadenwiderstand von maximal 
99999 Ohm, verianderlich in Schritten von 1 Ohm. R, kann ent- 
weder parallel zur Kapazitit C,(R? wie in Fig. 1) oder in Serie 
(R%) zu C, geschaltet werden. Als Nullindikator ist ein Kathoden- 
strahloszillograph verwendet worden. Die Wechselspannung wird 
der Briicke iiber einen symmetrischen Transformator zugefiihrt. Die 
am Eiskondensator hegende Spannung von 1 Volt wurde vor jeder 


Kapazitiitsmessung kontrolliert. 
x °e 
C, 





Fig. 1. 
Messbriicke. 


Die Messung der Kapazitit erfolgt durch Substitution. Zuerst 
wurde der leere Kristalltopf parallel zu C, geschaltet und der Wert 
C’, abgelesen, der fiir das Gleichgewicht der Briicke eingestellt 
werden musste. Dabei soll R?=oo oder also RS =0 sein. Dann 
wurde mit eingesetzten Eiskondensatoren gemessen und die Ein- 
stellungen C’, des Eichkondensators sowie der Wert R? oder RS 
abgelesen. . 

Unter der gemachten Voraussetzung C; = C, ergibt eine einfache 
Rechnung die Gleichungen fiir die Bestimmung von C, und R,. Im 
Falle der Parallelschaltung von R, und C, gilt: 
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und im Falle der Serieschaltung von R, und C, gilt: 


C+ Cz ” , 
C,=—4 4,5 —4-(4,--¢ 
7 1+0° 0? Re? i a 


sa Joe 
~ Cy+ C5 Ri (1 yi wo” Of RY ) 
Dabei bedeutet w = 2 zy die Kreisfrequenz und (C,—C,’) die Topf- 
kapazitit. 
Der Phasenwinkel wurde dann mit Hilfe der Gleichung 


1 
tg d= o RC, 
ausgerechnet. 

Der Flissigkeitsthermostat, den schon verschiedene Autoren 
beschrieben haben (vgl. z. B.8)), wurde jeweilen auf eine bestimmte 
Temperatur einreguliert und dann erfolgten die Messungen an den 
drei in den Kristalltopf eingesetzten Eiskristallen in Funktion der 
Frequenz. Die Bestimmung der Temperatur erfolgte mit einem 
Kupfer-Konstantan Thermoelement und einem Kompensations- 
apparat, wobei ein empfindliches Galvanometer als Nullinstrument 
diente. Die absolute Genauigkeit der Temperaturmessung betrug 
ca. + 0,5° C, die relative Genauigkeit mindestens 0,05° C. 


III. Die Kristalle. 


Die Ziichtungs- und Bearbeitungs-Methoden grosser Eiseinkri- 
stalle sind in einer friiheren Veréffentlichung ausfiihrlich beschrie- 
ben worden®). Die Zucht und die Verarbeitung konnten in einem 
Raum des Tiefkiihlhauses der Brauerei Hiirlimann AG., Ziirich, bei 
einer Temperatur von etwa — 20° C ausgefiihrt werden*). Aus den 
Einkristallen wurden quadratische Platten von etwa 2 cm Kanten- 
lange und 0,83—0,5 cm Dicke hergestellt. Jede einzelne Platte wurde 
mit Hilfe eines Polarisationsmikroskopes so orientiert, dass die 
Plattenebene entweder senkrecht oder dann parallel zur optischen 
c-Achse liegt. Die Genauigkeit der Orientierung schwankt, je nach 
Platte, zwischen + 1° und + 2°. 

Die schwierigste Operation besteht im Anbringen von Elektroden 
auf die orientierten, geschliffenen Eisplattchen. Da bekanntlich die 
DK-Messungen durch Fremdschichten (Klebstoffe u. dgl.) zwischen 
Dielektrikum und Elektrode sehr stark gefailscht werden, besonders 


*) Wir méchten Herrn Dr. H. HURL ANN an dieser Stelle unseren herzlichen 
Dank aussprechen. 
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wenn die DK des Dielektrikums gross ist, wie im Falle des Eises, 
war es nétig die Elektroden klebmittelfrei anzubringen. Das Auf- 
dampfen einer Metallschicht im Hochvakuum kam nicht in Frage. 
Es blieb daher nur die Méglichkeit eine Metallfolie anzufrieren. Eine 
diinne Silberfolie von ca. 0,4 « Dicke wurde auf einer ebenen Metall- 
platte auf 0°—1° C geheizt und das Eisplattchen daraufgelegt. Durch 
sorgfaltige Handhabung kann man erreichen, dass eine extrem 
diinne Schicht des Eises schmilzt, worauf die geheizte Metallplatte 
entfernt und das Eisplattchen wieder abgekihlt wird. Optische 
Untersuchungen im polarisierten Licht mit dem Anfrieren diinner 
Glasplattchen auf den Oberflichen von Eiskristallen zeigten®), dass 
die Regelation der geschmolzenen Schicht die Homogenitit des 
Eiskristalles nicht stért, so dass eventuelle Zwischenschichten von 
polykristallinem Eis nicht zu befiirchten waren. 
Die DK des Eisplittchens wurde aus der Beziehung 


re. 


a © and 
berechnet, wobei ¢« die DK, F die Plattchenoberflache und d die 
Plattchendicke bedeuten. 
Der maximale Fehler fiir F betrigt + 0,5%, fiir d etwa 2%, fiir 
C, etwa 1—2%. Ein verniinftiger Wert fiir die Genauigkeit der 


absoluten DK-Bestimmung ist deshalb von der Gréssenordnung 
3,5—4,5°%. Die Randeffekte wurden daher vernachlissigt. 


IV. Messergebnisse und Diskussion. 


Die Eiskristalle wurden zunichst aus gewéhnlichem Leitungs- 
wasser geziichtet, das keinerlei Reinigungsverfahren unterworfen 
war. Die Kristalle erschienen alle glasklar und optisch homogen. 
Die Beobachtung mit dem Mikroskop zeigte keine Lufteinschliisse 
und keine inneren Spalten. 


Bei den Messungen ergab sich, dass das dielektrische Verhalten 
von Kristall zu Kristall sehr verschieden war. Die gréssten DK- 
Werte, gemessen bei 50 Hz und — 5° C, schwanken zwischen 400 
und 900, der maximale Phasenwinkel zwischen 75° und 82°. Auch 
ist keine Regelmissigkeit in der Differenz zwischen der DK e,, 
welche im folgenden die gemessene DK bedeutet, bei welcher das 
elektrische Feld parallel zur c-Achse verliuft, und der DK e¢, , zu 
erkennen. (¢, , bedeutet analog die gemessene DK, wenn das elek- 
trische Feld senkrecht zur c-Achse verliuft.) Das Resultat einer 
Messreihe an zwei Plittchen aus verschiedenen Zuchten und ver- 
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schiedener Orientierung ist als Beispiel in Fig. 2 wiedergegeben. 
Im oberen Diagramm sind auf der Ordinate die gemessenen DK- 
Werte «, und ¢,,, auf der Abszisse die Frequenzen aufgetragen. 
Parameter der Kurvenschar ist die Temperatur. Im unteren Dia- 
gramm der Fig. 2 ist der Phasenwinkel 6, in Funktion der Frequenz 
und mit der Temperatur als Parameter dargestellt. 
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Fig. 2. 


DK und Verlustwinkel in Funktion der Frequenz. (Unreines Eis.) 
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Dieses Verhalten ist, nach den vorhin besprochenen Erfahrungen 
von GRANIER?) und von Smytx und Hitrcncock’), nicht schwer 
zu verstehen. Unsere Kristalle enthielten offenbar Verunreinigun- 
gen, die die elektrische Leitfaihigkeit und die DK bei tiefen Fre- 
quenzen sehr stark beeinflussen. Uber den genauen Mechanismus 

































































Fig. 3. 
DK und Verlustwinkel in Funktion der Frequenz. (Reineres Eis.) 


des Einflusses von Verunreinigungen auf das dielektrische Verhalten 
des Eises fehlen bisher noch die entsprechenden Experimente. 
Doch glauben wir nicht, dass Raumladungen irgendwelcher Art zur 
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Erklirung dieses Einflusses in unserem Falle herangezogen werden 
kénnen. So scheint uns z. B. unméglich, dass die Raumladungen 
dem wechselnden elektrischen Felde von 500 Hz noch folgen kénnen 
und doch ist die DK fiir diese Frequenz und die Temperatur — 10° C 
noch immer von der Gréssenordnung 200 (siehe Fig. 2). 

Wir versuchten danach, reinere Kristalle herzustellen. Das Wasser 
wurde zweimal destilliert, das erste Mal in Glas-, das zweite Mal in 
dicht abgeschlossenen Polyathylen-Gefissen*). Es war uns leider 


ELe-e-- 


PARAMETER 


FREQUENZ in KHe 








-50 -40 
T,°C 


Fig. 4. 
Thermischer Verlauf der DK. 


aus Zeitmangel nicht méglich, eine Leitfahigkeitsmessung des ge- 
reinigten Wassers auszufiihren. 

Die Resultate unserer Messungen an orientierten Eispliattchen, 
die aus diesem Wasser geziichtet wurden, sind in Fig. 3 graphisch 
wiedergegeben. Im oberen Diagramm sind die DK e, und «,, als 
Funktion der Frequenz mit der Temperatur als Parameter dar- 
gestellt. Es ist deutlich ein innerhalb der Fehlergrenzen reprodu- 


*) Die Destillation des Wassers verdanken wir Herrn Dr. H. GRANICHER. 
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zierbarer Unterschied zwischen ¢«, und ¢,, zu erkennen. Dieser 
Unterschied, d. h. die Anisotropie der DK des Eises, ist, wie auch 
kiirzlich von Powxss?°) vorausgesagt, sehr klein. Die relative Dif- 
ferenz (€.—€, ,)/€- fiir 1 kHz betragt etwa 15% bei —5° C, etwa 12% 
bei — 20° C und etwa 8% bei — 40° C. Der héchste gemessene Wert 
der DK, bei —5°C und 50 Hz, betragt 115 fiir e, und 100 fiir e, ,. 


Im unteren Diagramm ist nur der Phasenwinkel 6, , aufgezeichnet, 
um das Diagramm nicht zu uniibersichtlich zu gestalten. Das Ver- 
halten von 6, ist jedoch véllig analog zu demjenigen von 6, ,; nur 
ist der maximale Wert des Phasenwinkels stets etwas grésser fiir 
6, (69°—70°) als fiir 6, , (67°—68°). 

In Fig. 4 ist die Temperaturabhangigkeit der DK e, und «¢ , , wie- 
dergegeben. Parameter der Kurvenschar ist die Frequenz r in kHz. 
Interessant in diesem Diagramm ist das Verhalten der DK bei 
hohen Temperaturen und tiefen Frequenzen. Verfolgt man z. B. 
die Kurve fiir 50 Hz in Richtung abnehmender Temperatur, so 
sieht man, dass nach einer anfainglichen Abnahme bis ca. — 10° C, 
die DK einen Anstieg bis zu etwa — 30° C zeigt, um dann in der iib- 
lichen Weise zu fallen. Dieses Verhalten wurde nach unserer Kennt- 
nis in keiner der zahlreichen Veréffentlichungen iiber Eis erwahnt. 
Die anfaingliche Abnahme der DK zwischen — 5° C und — 10° C 
scheint nicht reell zu sein, sie ist auf Verunreinigungen zuriickzu- 
fiihren (vgl. auch das Cole-Halbkreisdiagramm Fig. 5 und 6). Hin- 
gegen scheint die Zunahme der DK von — 10° C bis etwa — 30° C 
reell zu sein. Betrachtet man niimlich die Temperaturabhiangigkeit 
der statischen DK ¢  (vgl. Fig. 10), so stellt man fest, dass &9 mit 
sinkender Temperatur zunimmt. Lasst man nun die Frequenz » 
vom Wert Null stetig zu héheren Werten tibergehen, so erscheint 
es einleuchtend, dass die DK fiir tiefe Frequenzen bei abnehmender 
Temperatur zuerst steigen und dann, infolge der Dispersion, fallen 
wird. 


Unsere DK-Werte fiir niedere Frequenzen und hohe Tempera- 
turen sind stets grésser als die DK des Wassers. Entgegen den 
Behauptungen einiger Autoren?)%)®), sehen wir keinen Grund, war- 
um die DK des Eises bei 0° C gleich der DK des Wassers sein sollte. 
Doch ist es méglich, dass unsere Kristalle immer noch Verunreini- 
gungen enthalten. Dies lasst sich beurteilen, wenn unsere Messungen 
in einem Cole-Halbkreisdiagramm!!) aufgetragen werden. 


Nach der Theorie von DreByz!?) wird die verallgemeinerte, kom- 
plexe DK durch die Gleichung 


e* =e’ — 18" (1) 
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dargestellt, wobei 


€" = (€)— Eco) - 


Hierin bedeuten: 


&9 die statische, ¢,, die bei sehr hohen Frequenzen (v > co) gemes- 
sene DK, w = 22» die Kreisfrequenz und t die Relaxationszeit. 
In der komplexen Ebene stellt Gleichung (1) einen Kreis dar, dessen 
Mittelpunkt auf der e«’-Achse liegt (e’ = (€ + €)/2) und dessen 
Radius gleich (&) — €,.)/2 ist. 

Unsere Briickenmessmethode liefert unmittelbar die relative DK, 
welche mit e’ von Debye identisch ist. e” findet man aus der Be- 
ziehung tg 6 = e”/e’. 

Die graphische Darstellung von e” = f(e’) ist fiir unsere Mes- 
sungen in den Figuren 5 bis 9 fiir fiinf verschiedene Temperaturen 








60 
e—- 


Fig. 5. 


wiedergegeben. Die Zahlen neben den Messpunkten geben die zu- 
gehérige Frequenz in kHz an. Fiir die Bestimmung des Halbkreis- 





Anisotropie der Dielektrizitatskonstante des Eises. 


T= -10,9 °C. 








T +-20,5 

















28 F. Humbel, F. Jona und P. Scherrer. 


mittelpunktes haben wir jeweils das Hauptgewicht auf die Mess- 
punkte bei hohen Frequenzen gelegt. Einen zusitzlichen Punkt zu 
diesem Zweck erhalten wir ferner aus ¢,,, wofiir wir den von LAMB 
und Turney!%) an polykristallinem Eis bei der Frequenz v = 
2,4-101° Hz (A = 1,25 cm) gemessenen Wert 3,18 oder auch den von 
Auty und Corn’) angegebenen Wert 3,1 zur Verfiigung haben. 


r 








C— 


Fig. 9. 


Die Diagramme fiir 7’ = — 5° C und T = — 10° C zeigen den ty- 
pischen Einfluss von Verunreinigungen bei tiefen Frequenzen: Fiir 
y kleiner als 1 kHz bzw. 300 Hz steigen die Verluste (e”) rapid an 
und die entsprechenden «’-Werte erscheinen auch zu gross’). Mit 
sinkender Temperatur nimmt auch die Frequenz ab, bei welcher 
die Messpunkte vom Halbkreis abzuweichen beginnen. Bei — 5° C 
ist diese Frequenz etwa 1 kHz; bei —10,9°C etwa 300 Hz. Bei 
— 20,5° C liegt der Messpunkt fiir die tiefste gemessene Frequenz 
(50 Hz) noch sehr gut auf dem Halbkreis. 

Der Schnittpunkt des Halbkreises mit der e’-Achse auf der Seite 
der tiefen Frequenzen ergibt die statische DK ¢). Man kann somit 
durch diese Extrapolation eine Auskunft tiber das thermische Ver- 
halten von ¢, bekommen. 

Fig. 10 gibt den Verlauf der statischen DK fiir die Richtungen 
parallel (¢,, leere Kreise) und senkrecht (e€,,,, volle Kreise) zur 
c-Achse wieder. Die Messpunkte fiir —5°C weichen von dem zu 
erwartenden Verlauf der Kurven ab, was méglicherweise einem 
Messfehler zuzuschreiben ist. Die Extrapolation der Messreihen er- 
gibt fiir die Temperatur 0° C die Werte e,, = 104,8 und ¢, ,, = 91,2. 
Nach Auty und Cots’) ist die statische DK des polykristallinen 
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reinen Eises 91,5 bei — 0,1° C. Es ist méglich, dass das von Auty 
und Coz untersuchte Eis einen grésseren Reinheitsgrad besass als 
unsere Kristalle, da diese Autoren diesem Punkte spezielle Auf- 
merksamkeit schenkten. In diesem Falle kénnte der Wert 91,5 als 
ein Mittelwert von ¢,, und ¢,,, angesehen und unsere Werte ent- 
sprechend auf das ,,reine Eis korrigiert werden. Es scheint uns 
aber andererseits nicht sicher genug festzustehen, ob das von Auty 
und Coe untersuchte Eis wirklich eine statistische Verteilung der 
Kristallitorientierung aufwies. Unsere Erfahrungen in der Zucht 
von Eiskristallen lehrten naimlich, dass die einzelnen Kristallite nicht 
klein sind (so dass ein Eiszylinder von etwa 8 cm Durchmesser und 
4 cm Hohe deren nur 6 bis 7 enthalten mag), und dass ferner zwei 
Eisblécke, die unter denselben Bedingungen geziichtet werden, 
meistens nicht die gleiche Einkristallkonfiguration aufweisen. AuTY 
und CoLs untersuchten aber die statistische Unordnung der Kristall- 
individuen nicht an demselben Eisstiick, das fiir die dielektrischen 
Messungen bestimmt war, sondern nur an einem anderen Stiick, 
das unter ahnlichen Bedingungen geziichtet worden war. Das ther- 
mische Verhalten der statischen DK gemessen von Auty und CoLE 
ist in Fig. 10 durch die Kreuze wiedergegeben und zeigt einen inner- 
halb der Fehlergrenzen ahnlichen Verlauf wie unsere Kurven. 


In Fig. 10 sind zum Vergleich auch noch die theoretisch von 
Pow tess?) ausgerechneten Werte der statischen DK des Eises auf- 
getragen. Die Rechnungen von PowtxEs machen von der allgemeinen 
Fréhlichschen Theorie!4) der statischen DK Gebrauch und stiitzen 
sich auf das Paulingsche Modell fiir die Struktur des Eises!5). Nach 
diesem Modell kann eine Gruppe von H,O-Molekiilen in einem Eis- 
kristall nur eine bestimmte Anzahl von Konfigurationen annehmen, 
die aus der Wechselwirkung zwischen den Molekiilen resultieren. 
In einem ersten Fall wird von PowLEs angenommen, wie PAULING 
vorgeschlagen hat, dass alle Konfigurationen die gleiche Wahr- 
scheinlichkeit besitzen. Dann ist die DK exakt auszurechnen, €&9 
wird 77,3 bei 0° C und man kann sogar die Anisotropie der DK als 
nicht gross“ voraussagen. Es ist aber nicht méglich dariiber eine 
‘quantitative Angabe zu machen. Die aus der Powles’schen Formel 
ausgerechneten Werte fiir diesen Fall sind in Fig. 10 mit den Kreuz- 
kreisen eingetragen (Powxss I). In einem zweiten Fall wird von 
Pow.es der Einfluss der elektrostatischen Wechselwirkung zwi- 
schen den Molekiilen auf die Wahrscheinlichkeit der Konfiguratio- 
nen mitberiicksichtigt. Die Genauigkeit des Resultates ist dann 
schwierig anzugeben und eine Aussage iiber die Anisotropie der 
DK ist nicht méglich. ¢) wird gleich 103 bei 0° C. Das thermische 
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Verhalten der statischen DK fiir diesen Fall ist in Fig. 10 mit den 
Dreiecken angegeben (Powtss II). 

Der Verlauf der Kurven unserer Messungen steht mit dem theo- 
retisch berechneten Kurvenverlauf von Powxs in guter Uberein- 
stimmung. 

Zum Schluss soll noch die Temperaturabhingigkeit der Relaxa- 
tionszeit tT aus unseren Daten berechnet werden. 


7 


Powles I 
Powles IT 


Auty - Cole 








“a a)! lUe.hlUle Cl 
7% 
Fig. 10. 
Thermischer Verlauf der statischen DK. 


Die Grésse t wurde erstmals von DEBYE?*) aus den Messergeb- 
nissen von Errera?) ausgerechnet. Aus den Gleichungen (2) und 
(3) ist leicht einzusehen, dass wenn 

& = (9 + &)/2 (4) 
d. h. wenn die DK gleich dem arithmetischen Mittel von ¢) und ¢, 


oe wt = (60 + 2)/(e + 2) (6) 
Die Kreisfrequenz , fiir welche Gleichung (4) gilt, kann dem Dia- 
gramm der Frequenzabhingigkeit der DK entnommen werden 
(Fig. 3). Somit wird t bei bekannten ¢9 und «,, mit Gleichung (5) 
ausgerechnet. Desyer findet fiir die Temperatur —5°C, t= 
2,7-10-* sec. Unsere Messungen ergeben, fiir dieselbe Temperatur, 
die Werte t, = 1,4:10-* sec und t, , = 1,7-10- sec. 
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Auty und Cor’) geben eine etwas abweichende Definition der 
Grésse t. Gemiiss den Gleichungen (1), (2) und (3) wird ein ,,Dipol- 
verlustfaktor“ tg gp als 

tgp =e" (e’ —e,) = or 
definiert. Traégt man tg y als Funktion der Frequenz auf, so erhalt 
man eine Gerade. Die Frequenz »,, fiir welche tg » = 1, ergibt dann 
die Relaxationszeit t aus der Beziehung t = 1/2 2 »,. Diese Relaxa- 
tionszeit unterscheidet sich von der obigen, nach DeBYE ausgerech- 
neten, um den Faktor (€) + 2)/(e.. + 2). Aury und Coxe erhalten 


=2 
ae 


BERECHNET 
NACH AUTY-COLE 


AUTY-COLE 


BERECHNET 
NACH DEBYE 


DEBYE -ERRERA 





o® T° 
36 37 38 39 40 4 42 43 46 45 








0-0 -2 -30 -40 -50 -60 =? T,°C 
Fig. 11. 
Temperaturabhangigkeit der Relaxationszeit. 


fiir die Temperatur —10,8° C den Wert t = 6,0-10-5 sec. Unsere 
Messungen ergeben fiir die gleiche Temperatur: 


t, = 5,8-10-5 sec tT, ¢ = 6,1-10- sec. 


Fig. 11 zeigt das Verhalten von log t in Funktion der reziproken 
Temperatur. Die Werte von Auty und Cots erfiillen die Beziehung 
t = A-eB/RT, Unsere Werte sind angegeben durch leere bzw. volle 
Kreise fiir die Relaxationszeiten in Richtung parallel und senkrecht 
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zur c-Achse berechnet nach DresByeg, ferner durch leere bzw. volle 
Dreiecke fiir die nach Aury und Cote berechneten Relaxations- 
zeiten. Fiir die Temperaturen unterhalb — 35°C wird die obige 
exponentielle Temperaturabhingigkeit von unseren Werten nicht 
mehr befolgt. Die Genauigkeit unserer iiber so viele Umrechnungen 
bestimmten t-Werte und das untersuchte Temperaturintervall sind 
nicht gross genug, um eine befriedigende Deutung dieses Verhaltens 
zu erméglichen, doch kénnte eventuell eine Temperaturabhangigkeit 
der Aktivierungsenergie B herangezogen werden, welche wahr- 
scheinlich von den Verunreinigungen beeinflusst wird. 

Dieses Verhalten der Relaxationszeit, ebenso wie die Temperatur- 
abhingigkeit der statischen DK, die im untersuchten Temperatur- 
intervall mit abnehmender Temperatur immer mehr ansteigt, ver- 
langt weitere eingehende und systematische Untersuchungen bis 
zu tieferen Temperaturen. Auch scheint uns leider der Einfluss von 
Verunreinigungen, trotz den letzten wertvollen Untersuchungen von 
Auty und Cots, noch nicht klar genug zu sein. Die vorliufig beste 
Erklarung fiir das eigenartige dielektrische Verhalten des Eises 
beruht ja bekanntlich auf dem Bild von molekularen Reorientie- 
rungen als Folge von Gitterdefekten!*)!”), die sich, wegen der 
Wasserstoffspriinge, leicht im Gitter bewegen kénnen. Die Verun- 
reinigungen erhéhen die Dichte dieser Gitterdefekte und somit die 
Polarisation. Die Frage, welches der tatsichliche Mechanismus und 
was der Einfluss verschiedenartiger Verunreinigungen sei, bleibt 
vorlaiufig noch unbeantwortet. 
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On the Divergence of Perturbation Theory for Quantized Fields 
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(12. X. 1952.) 


Summary. The convergence of perturbation theory is investigated for a scalar 
field with an interaction term A y*. It is shown that for a certain energy region 
the series diverges for all values of /. 


1. Introduction. 


The most important progress in field theory in the last years was 
that it was possible for certain theories to absorb the infinities oc- 
curing in the perturbation theory in renormalization constants. The 
possibility of making a theory finite by renormalization stems from 
the fact that no new types of infinities appear as one proceeds to 
higher approximations. Thus the infinities can be absorbed by a 
finite number of compensating terms in the Lagrangian. 


A renormalizable theory would provide a complete mathematical 
scheme, if the power series in the coupling constant of perturbation 
theory converged after renormalization. On the other hand, if the 
series proves to be divergent, there seems to be no significance in the 
fact that higher orders do not produce new infinities. 


Until some time ago it was the general believe that at least for 
small values for the coupling constant the series converges. This 
would suffice, as this function of the coupling constant could be 
extended by analytic continuation to any desired value for the 
coupling constant. Just recently Dyson found a physical argument 
which makes it plausible that there is no radius of convergence for 
perturbation theory in quantum electrodynamics. Yet no mathe- 
matical analysis of this problem has been given so far. The expres- 
sions of higher order terms in perturbation theory become so complic- 
ated that nobody has been able to show whether terms of arbitrary 
high order are actually small. The purpose of the present paper is 
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to study this problem for a simplified model!). It seems reasonable 
to restrict oneself to scalar fields as it does not seem likely that 
the question of convergence is tied up with the transformation pro- 
perty of the field?). 

This leads him to consider two scalar fields ® and y which are 
coupled by an interaction term 4 ® y?, being linear in one field 
and bilinear in the other. A closer investigation of this theory shows 
that as far as convergence is concerned this theory behaves like a 
scalar field with a non-linear term / y?. 

It can be seen easily that this leads to a renormalizable theory 
containing a renormalization of mass and field strength. The result 
of the following analysis will be that the series in powers of A for 
the propagation function of the nonlinear field diverges within the 
energy region where no new particles can be created. Thus it is 
essential that there are no particles with vanishing mass. The 
special form of the interaction term does not seem to be very im- 
portant as long as it is hermitian’). 

The proof of the divergence rests essentially on the fact that cer- 
tain terms always have the same sign as long as the interaction term 
is hermitian. The same property is used in Schwinger’s proof that 
the charge renormalization always decreases the charge. This seems 
to indicate that the divergence of pertubation theory is a deep 
lying feature of field theories which cannot be removed by mixing 
several fields. 

The question whether solutions exist other than expansions in 
powers of the coupling constant will not be considered in this paper. 
Preliminary investigations on this point seem to indicate that not 


1) It was observed independently by Dr.C. A. Hurst and by the author that it 
is possible to determine the convergence behaviour of perturbation theory for such 
a model. The work of Hurst is contained in his (unpublished) thesis Cambridge, 
January 1952, a paper of his is to appear in the October 1952 issue of the Proc. 
Camb. Soc. The author announced his results at the Copenhagen conference, 
June 1952. It turned out that our results parallelled each other, though both works 
were incomplete in some respect. A mistake in the (unpublished) work of the author 
was pointed out to him by Dr. G. KAtiin. The work of Hurst does not deal with 
renormalization and contains a mistake in sign (in formular 23). But the short- 
comings of our works did not overlap and in view of the importance of the subject 
the author thought it worth to give a complete proof. The present analysis uses 
in one point an idea of Hurst and is a simplification over both original works. 

*) One may be inclined to simplify the problem by reducing the number of the 
four dimensions and investigate a one-dimensional problem. But for this case one 
can easily construct a convergent theory. The divergent momentum integrals seem 
to be connected with the 4-dimensional structure of space-time. 

3) One could carry out the following investigations along the same lines with 
an interaction term 4 y* though this would be somewhat more involved. 
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only a power series in powers of A but a power series in powers of 
4 — dy diverges as well. This would mean that if there are solutions 
at all they must be non analytic for all real values of 4. There exists, 
of course, the possibility that for certain values for 2 solutions do 
exist even for the unrenomalized equations. But it seems very un- 
likely that for e?/4 2 = 1/187 there exists a solution for quantum 
electrodynamics. Particles of spin 4 and spin 0 have the same 
charge and it is hard to believe that these two different theories 
have solutions for the very same coupling constant. 


2. Renormalization of the theory. 


At first we shall briefly outline our notation. We use real world 
coordinates with the metric go = 1, 911 = Joo = J33 = —1. The 
vector product of 4 vectors shall be written in the simple way pk = 
Poky—P1i ki — Pok, — pgks. The operator P effects the chronological 
ordering of the following factors. We further use Wicks symbols for 
the ordered product: ...: where the expression within the double 
columns has to be ordered in such a way that the positive frequency 
parts of all operators stand to the right hand side of the negative 
frequency parts. For the Dalembertian we write 0? = 0;0%. For 
the other symbols we use the standardised notations. 

We shall now consider a scalar theory characterised by the La- 
grangian 


L(x) = 5:0, p 0! y—m? y?: (2,1) 


from which emerge the field equations and the commutation rela- 


tions 
(O?+m?)y=0 = [yp(z), p(2’)]}=—tA(a—z’). — (2,2) 


The only other expression we need for the following investigation is 
the vacuum expectation value for chronologically ordered operators 


P(0| y(a) v(2’)|0) = —i f ap (2.8) 
Now we introduce an interaction term 
L"* (2) = A: p(x): (2,4)?) 


which we shall treat as a perturbation in an interaction representa- 


ep (4—2’) 


m? — p * 


1) Note that: y*: is hermitian. If one takes y*® as interaction term one gets 
spurious renormalizations owing to the nonvanishing vacuum expectation value 
of the interaction term. One would obtain graphs of the type |—o which give an 
additional mass renormalization. 
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tion. Ais the coupling constant which has a dimension of a reciprocal 
length. On calculating the S-matrix with the interaction term (1,4) 


S'=14 id [ dey: y(l):4+ 3, (14)? P [dary dy: y8(1):: y?(2): te 


1 a 
A m= : : (2,5) 


one gets graphs where every vertex carries three lines quite similar 
to quantum electrodynamics. The differences between the graphs 
in this theory and quantum electrodynamics are: 

a) some graphs are identical in this theory while graphs of the 
same topological structure are different in quantum electrodyna- 
mics because we have just one type of field. 

b) For all internal lines we have to insert the propagation factor 
—1/m? — p?. 

c) There are no factors steming from the vertices nor are there 
any trace calculations. 

One can now easily investigate in which graphs diverging integrals 
occur over energy momentum vectors of virtual particles. Follow- 
ing the analysis which Dyson has carried out in quantum electro- 
dynamics we can find the primitive divergent graphs simply by 
comparing the powers of the energy momentum vectors in the nu- 
merator and denominator of the expression for the graph. Because 
from each vertex there emerge three lines, the number of external 
lines |, and the number of internal lines /; are related by 


21,+1,=8n. (2,6) 


Each internal line contributes a denominator with the squared 
energy momentum vector to the graph. The number of energy mo- 
mentum vectors over which one has to integrate is equal to 1; —n+1. 
The difference K between the powers of denominator and numerator 
is therefore K = 21,—4(l,—n+ 1). If we use (2,6) this becomes 
K=n+1,—4. In the following table we put K in dependence 
of n and lI,. 





| a 
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K is less than 1 only for those values of n and |, indicated by the 
heavy line. The graphs with |, = 0 sum up to a common phase factor 
giving the probability that vacuum remains vacuum. They shall 
not be considered further. The vacuum expectation value of y in 
the HrisenBERG representation does not vanish in this theory and, 
therefore, these are non-vanishing graphs with just one external 
line. The graph n = 21, =1 does not exist but only the graph 
n=31,=1--— >. The last diverging graph is the graph |, = 2 
n = 2——C>— which corresponds to the selfenergy graph in quan- 
tum electrodynamics. 


It seems to be interesting to consider a regularized theory where 
we replace the propagation factor by 1/(m?— p*)*. In this case 
K = 4 (n—1) and, therefore, there are no divergent graphs. Yet it 
turns out that even for a propagation factor 1/(m?—p?)*, t being an 
arbitrary integer, the series of perturbation theory diverges. There- 
fore, in a theory with a form factor which changes the propagation 
function only in a mild way of the above form, there is but little 
hope for convergence of perturbation theory. 


Our theory differs from quantum electrodynamics in as far as n 
does not disappear in the expression for K. In our theory types of 
graphs which are divergent in lower orders become convergent in 
higher orders?). 

Though these graphs converge they contain finite renormalizations 
which have to be separated out as well. 


We determine the renormalization constants by the condition that 
higher order corrections vanish in the non-relativistic limit. This is 
fulfilled if we cast the propagation function in the general form 


fo — [1+ C(m?—p?)] with C(o)=0. (2,7) 


m?— 


We next consider sub-graphs which are only connected by one 
line with the main graph. As we have seen the lowest order graph 
of this type is divergent. It turns out that inserting these sub- 
graphs in any line of a graph amounts merely to a mass renormaliza- 
tion. On account of energy momentum conservation the line joining 


1) This seems to be connected with the fact that the coupling constant has the 
dimension of a reciprocal length. In renormalizable theories where the coupling 
constant is dimensionless, n cancels out in the expression for K. Theories where 
the coupling constant has the dimension of a length are not renormalizable and 
in this case n remains with the negative sign in K. Thus our theory may be called 
super-renormalizable. 
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the sub-graph and the main graph has zero energy momentum 
vector. One can see this by the following simple example: 





Fig. 1. 


The effect of inserting this sub-graph is doubling the propagation 
factor — 1/m? — p? and multiplying the whole expression by 


-i f ; 7 -i 
D3= =, / dk dk ( i) mk? m= (k+ hE’ 
The effect of this insertion is exactly compensated by adding a term 
Dg: y?: to the LAGRANGIAN. 

There is no finite correction left over after renormalizing this 
graph. The same holds, of course, if the sub-graph is of a more 
complicated structure. A closer consideration shows that all graphs 
which contain sub-graphs, joined to the main graph by only one 
line, are exactly cancelled if one adds a term D: y?: to the La- 
GRANGIAN. D is the sum of all graphs containing one external line 
only. Henceforth we shall disregard graphs containing this type of 
sub-graphs. 

The renormalization of self-energy graphs can be carried out in the 
same way as in quantum electrodynamics. We call a self-energy part 
proper, if it cannot be divided into two parts by cutting one line. 
Then the propagation function 4, for the non-linear field which 
contains the effects of all self-energy graphs can be expressed by 


A, = A, + 4, ZA, (2,8) 
_ where 2 is the contribution from all proper self-energy graphs. In 


momentum space 2' depends only on the invariant p? and we ex- 
pand 2 in powers p* — m2. Then the above expression becomes 


Ap(p?) = Ap(p?) + Ap(p?) [2 (m?) — (m? — p?) 2” (m?) + 
+ (m?—p?)? R(p?)| Ap (p?) . 
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The first term can be compensated by a mass renormalization‘) term 
in the Lacrancian 2(m?):y?:. The second term has to be com- 
pensated by an addition 


=’ (m2): 0; yt yp —m? y?: (2.9) 


to the Lacrancian. It represents a renormalization of the field 
strength”). After separating out these two terms the remaining 
expression A, is of the desired form (2,7). 

Renormalization works here in a similar though in a somewhat 
simpler way as in quantum electrodynamics. In the next section 
we shall investigate the expression for X more closely. Perturbation 
theory gives an expansion for 2 (p?) and our aim will be to show 
that the series diverges for a certain domain of p?. 


3. General expressions for A;. 


We now consider the contribution of a self-energy graph of n™ 


order to 44 ——: 


We forget at the moment about renormalization. Then the contri- 

bution will be ; : 
—% a) 

m?— p? ~n (p*) m3 — p? * 


+ 


2, consists of 3 n/2—1 internal lines each contributing a factor 
— i/m? — k?, if k is the energy momentum vector of the correspond- 
ing line. There are n/2 energy momentum vectors over which one 
has to integrate, the others are fixed by conservation laws. For 
carrying out the integration over n/2 k’s we combine the denomi- 
nators in the expression for 2, with the identity 


Tn-1 


1 XY 
1 = (,47;-1)! * ' 
sap ae Mean JO ee 
sli - 0 0 Gj 


0 
(1 ay)"*~* (ay — y)""71 (aq)! 


[Ay + % (Ay—p)+ +++ +2 (An —4y_y)]**7% 





(3,1) 


We shall use the abbreviated notation [dz for the integral over 


1) The compensating terms given here are not quite correct. This will not be 
analysed, however, as it has no influence on the results of this paper. 

*) There is a formal complication as this term in the Lagrangian yields normal- 
dependent terms in the Hamiltonian. This point will not be discussed here as it 
has no bearing on our main problem. 
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the x... 2, including the normalization factor n! and the (1 —2,)”... 
(x,)""-1. This casts the expression for 2 into the form 


oy3nj2—1 [ * dk, dk n/2 
(—1) a ae gee 
x [m?— (k:Q (5; +217; (x) pk; +s (x) p*)y* m[2+1 (3,2) 


The quadratic form Q;;, 7; and s depend on the « only. Because 
m? occurs in each denominator, it is not multiplied by an x-depen- 
dent factor. For carrying out the integration over k, we have to get 
rid of the terms with (pk,). This is accomplished by the substitution 


° - —1 » 
kj >k;—Q3 rep. 


Here we consider @ as a matrix of the rank n/2 and Q~-1 is the reci- 
procal matrix. This transforms the bracket of the integral into 

( ) —_> kh, Qik; cL 2(s — "Qi; r;). 
Because Q is a symmetric quadratic form with real coeficients it 
can be made diagonal by a transformation with a determinant 1. 
After this transformation we are left with an integral 


sii ‘fae fa ae — [m?—B;(x) k?—p?a(x)]-"?** (3,3) 


where f; are the Eigenvalues of Q and « (#) = s(x) — r,(x) Q;;(a)r,(a). 
The integral over k will be carried out with the aid of the well-known 
formula 

[ dk im? 


[L+ak}+2 ~ aL n(n+1) * 
This leads us to the result 


; vane _n (n/2-2)! / da 1 
yl (p?) aa (—1) . (4 7) (3 n/2—2)! / |Q(x)?| [m2 — p? a (x) ]”/2-2 , 


where Q is the product of the Eigenvalues of Q,;. The integral, there- 
fore, does not depend on the Eigenvalues separately, but on the 
determinant only. The Eigenvalues of Q can be shown to be positive. 
(Appendix a.) An important property of «is o < a(a#) <1. 


This comes from the fact that there are no real processes for 
p? < m*. The creation of new particles by an external field with 
the energy momentum p is given by the real part of 2,(p?). The 
integral for 2’ is real as long as there are no zeros of the denomi- 
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nator. If poles occur in the integrand one has to take the path of 
integration in complience with the path of integration for the 4,- 
function. It amounts to adding a small negative imaginary part to 
the mass and taking the principle values for all integrations. As no 
particles can be created for p? < m? there cannot be a real part of 
+ and, therefore, no poles in the integrand. A direct proof will be 
given in appendix a. 

A graph of n'" order of the form —- 2, —— provides us with 
a function of p? which is of the following form: 


—i (42)-" — (n/ =O * dx 1 ; 
mip [meat me) [Ol TwPprapee | (OM) 
Here we have taken into account the factor 1” of the general expres- 
sion (2,5). We note that the factor in the square bracket is positive 
definite for p? < m?. One can see from general considerations?) 
that —1 & (p?)/m? — p? is positive definite, so that the sum of the 
contributions of all graphs is bound to be positive (within this 
energy region). Our result tells us that the contribution of each 
graph separately is positive. 

We now consider an irreducible self-energy graph and carry out 
renormalization according to the prescription given in the preceding 
section. The expression given by an irreducible self-energy graph 
will be of the general form (3,4). We first separate out the mass re- 
normalization. 

+ (et. 3 1 .. 
] [Q|? | (m= p? a)/2-2 [m2 (1—a) 2-2 J 


m2 — p> | 


n \ ” dz 1 se 
=— e—))/ dy, / |Q|? [m2 (1 — a) + (m?— p?) ya]? * 3,5) 
0 


In the case of n = 2 one has to renormalize before carrying out the 
k-integration. The result is again of the form (3,5). A second sub- 
traction allowing for the renormalization of the field strength gives us 


aie { 1 . Res: eer, 
(- 1) an | [m2 (1—a) + (m?—p®) aysJ"/® — [m® (1a)? J 


1 - 


(m?—p*) y ‘ 
ba (= “Hy 2 | dun dy yo |G ic a” [m2 i-wsee enae : (3,6) 


1) KALLEN, Helv. Phys. Acta, 25, 417 (1952). 
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After renormalization the contribution of an irreducible self-energy 





i m? — n/2! dx dy a? y, 
m2 — p? i ae (3n/2—2)! | aria a | - (8,7) 


After the double renormalization the expression in the square 
bracket is again positive for p? < m?. 


If one deals with reducible graphs, which are graphs containing 
subgraphs of the self-energy type, renormalization replaces the 
subgraph by a propagation function of the form (3,7). Thus one has 
to insert into the skeleton graph instead of the propagation function 
— 1/m? — p? a more general function of the form 


1 
< >b>0. 


; aig ‘ 3 
/ Iman peppA with 1>a> * 


Inserting this propagation function into the lines of an irreducible 
self-energy graph leads us back to an expression of the type (3, 4). 


There is a difficulty arising, if we encounter self-energy graphs 
of the form —o—o—, where two or more self-energy graphs are one 
on top of the other. In this case one cannot use the line on the right 
side and on the left side of a self-energy graph in the process of 
renormalization and there remains m*— p? in the numerator. Then 
one has to apply more general integration formulas, a case which 
will be studied in the appendix b. This difficulty would not arise, 
if one considered a regularized theory with a propagation factor 
— i|(m? — p®)2. 

The result of this section is that the expression for A,, is of the 
form — i[1+2G(p?)|/m* — p? where the sum runs over all self- 
energy graphs which are renormalized according to the general 
pattern. The contribution of each graph is positive for p? < m?. 
This gives us the advantage that we get a lower bound for G by 
picking out certain terms of the sum which are easy to calculate. 
This will be done in the next section where we give an estimation 
of certain irreducible self-energy graphs. Our aim is to find self- 
energy graphs of n order which are of such a simple structure that 
they can be well estimated, yet their number should be large enough 
to effect the divergence of the series. 
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4, Estimation of terms of n‘" order. 


For estimating contributions of n™ order to G we have to look 
for graphs for which one can find a lower bound for «, which is 
greater than nothing. This leads more or less uniquely to considering 
graphs of the following type: 








Fig. 2. 


The part within the broken lines is a vertex part K of the order 
n — 8. Of these graphs we take those in which the line with the 
energy-momentum-vector k enters in the middle of the graph and 
where the n/2 — 2 points on the left hand side are joined with the 
n/2 — 2 points on the right hand side: 














Fig. 3. 


The number of those vertex parts is equal to (n/2 — 2)! corre- 
sponding to the possibilities of joining the n/2— 2 points on the 
left hand side with the points on the right hand side’). Here and 
in the following we consider only topologically different graphs as 
the number of connected graphs which differ only by the labelling 
of the points equals n!. This cancels exactly the n! in the denomi- 
nator of our general expression. The total number of contributions 
of the n order is less than 


e 7 7 3 ! 
(3n—1) (8n—38)....1 cama ; (4,1) 


1) Here we disregard the possibility of choosing the three y’s of each vertex 
in different ways. This gives us a factor of the order 3°"/2 which is, however, of 
no importance in the question of convergence. 
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This is the number of possible pairings of 3 y-factors. The: 
number of graphs of n‘ order is less than that, as pairings of y’* 
belonging to the same point do not contribute. This number is to 
be compared with the number of graphs we have selected. This is 
equal to n! (n/2—2)!, being about the 8-3"? 2-"? th part of the 
total number. For calculating K we use the formula (3,1). The de- 
nominator will be of the general form 


m?—k,Q;;k;—2k, (ai r+a,k)—b,r?—2 bork—bgk?. (4,2) 
Q and a and b depending on « only. After carrying out the integra- 


tion over dk,....dknjo-5 we arrive at the following expression: 


ae ye (n/2—3)! [4, 
K =: (4 zt) +4 ean | dx x 


x [m2—d, 12 —de r k— dg k2}-"!2+2 | Q (a) |-2 (4,3), 


d, = b, — ai Q;;° a 

da = by — a, Qj; a (4,4) 

ds = bs — a3, Q;;' a}. 
The condition that there are no real processes for k = 0, r? < m* 
and for r = 0, k2 < m? and r = —k, k? < 0 supplies us the inequa- 
lities 

d,+d, ss 
0<d,<1 <-_<' ee (4,5) 


If one combines this expression with the external line carrying the 
energy momentum vector k and carrying out the integration over 
k one is led by a similar argument to the inequality 


2 < dy dg. (4,6) 


When combining this expression with the denominators for the re- 
maining lines of the complete graph fig. 5 we first take those with 
the energy momentum vector r and r + k: 


. es 
m*—r2 m*—(r+k) 


a —n (n/2—3)! 
» K= (4a) Spa 


[ de dy x 


ase. eee 
4 ~ [m2 — (c, r2 +2 cy rk+ cy k?)]"/? 
with 


Cy=1—yo(l—d;)  ¢2=Yi—Yo(1—de) C3 = Yi — Yo (1 —ds) 
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In view of the inequalities (4,6) we get 


0<e<1 c, << 45% 
0<¢,< 1 O< G@< qty. 


In (4,7) we have used the abbreviated notation / dy for 


1 nn 


2) (5 a 1) [ay, / dys yt?-8, 


0 0 


We join the remaining denominators to the expression (4,7) with 
the help of three auxiliary variables z. 


[m—(p+k)2}-1 [m2(p+r+k)?2]-3 [m2k?]-1 [m?—c, 1-2 cr kc k2]}-"? 
sae Nan See (4,10) 


Then we can carry out the integration over r and k and arrive at 
an expression 


in 
> 3}! ‘ -21¢Q’|_— 
(3 ) | dxdydz a) —— 


(42)-" (4,11) 


=e _9)n(=+2)(5-+3)- = (2-1) [m?—ap"] 5-1 


The quadratic form of the denominators before integrating over r 
and k is given by 





k 


1-2, 


1—2,+ 2C, 














From this we calculate the expressions 


Q’ = (1—2g + 23 €3) (21 —2_ + 231) — (21 — 22 + 2g Ce)? 
a={(1—2,) [(1— 2g + 2g €3) (21 —2_ + 23 C1) — (41 —2q + 2 Co)?] — 
— (1 — 2g)? (2) — 2g + 23 1) + 2(1—2q) (@1—2e) (21 —22— 23 C3) — 
— (2) — 2g)? (1— zg + 2g) } { (1 —2g + 2 €g) (2 — 2% + 23 ¢1) — 

— (2;—2q + 2 Ce)? }-?. (4,12) 


The last step is carrying out the double renormalization and multi- 
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plying with the last denominator which contains the energy momen- 
tum vector p only. Thus we are finally led to 


[ _dadydzdu 
|Q (x) |?/Q’|2 


yy a2 (m* — p*) 
[m? (L—a)+ 0 ay ly (me ph * 





(4,13) 


In this expression the y- and z-dependence of « and Q’ is known 
whilst the z-dependence is unknown. The minimum-value of a, 
however, with respect to the 3n/2 —6 variables x is greater than 
zero in view of (4,9). We obtain a suitable lower bound for H by 
replacing all z-dependent quantities except Q in the numerator and 
denominator by their minimum or maximum value respectively. 
Then we can break the integral for H into two parts: 





H >(3/ (n/2 — 3)! (n/2)! / dz dy du u, «? (m?— p*)min 
4a} (3n/2—7)! (n/2+2)! | Q’ |? {m2 (1—a) +0 w, Up (m2— p2)}"/2+1 


x / da |Q(a)|-2. (4,14) 


Furtheron we concentrate on estimating the integral over dx. This 
can be done by giving a close lower bound for (4,3). For this goal 
we observe that in the region p? < m? there are no displaced poles 
(in the terminology of Dyson, Phys. Rev. 75, 1736 (1949) in the 
ko-integrations. Thus we put in (4,2) r = k = 0 and rotate the path 
of integration for each kp-variable to the imaginary axis. Then each 
propagation factor becomes positive definite and we can use the 
fact that the geometric mean value is less than the arithmetic mean 
value for combining the denominators?) for K: 


1 1 (3 )3* 
m>— ki m—ks, (3tm?—ki—---—Ky,)** - 








(4,15) 


Here and in the following we put n/2 —2 = t. We obtain a lower 
bound for (4,3) by using (4,15) instead of (8,1). It proves useful not 
to use the conservation laws from the vertices for expressing 2 t 


1) This method is essentially due to Hurst: (Compare introduction). 
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k-variables as linear combinations of the basic t k-variables. We 
postpone this and write the k-Integral with matrix notation 


. (4.28 
(34)%* / d430K Sint Kn . (4,16) 
The 4-3t-dimensional vector K is composed oi the 3¢ fourvectors k;. 
The 4-2¢-fold 6-function expresses the conservation-laws stemming 
from the 2t vertices. They are written in matrix notation A,,k; = 0, 
where A is a 4-2t by 4°3¢ rectangular matrix. Each row of A con- 
tains at most three times + 1 and 0 everywhere else. This is be- 
cause there are three lines at each vertex and the conservation law 
puts the sum of the three energy-momentum-vectors equal to 0. 
The 4-2¢ rows of A may be considered as a system of 4-2¢ linear 
independent vectors a; which span a subspace S. The 6-function 
demands that K must be orthogonal to all a, or that the integration 
has to be carried out over the complementary subspace T of 4-¢ di- 
mensions. K can be split into a component Ky, in S and a component 
Ky, in T. For the domain of integration K? = Kj holds. In order 
to eliminate the 6-functions we choose an orthogonal set of 4-t vec- 
tors n; in T' and make the linear transformation 

4-3¢ 

N 4-2t 


MK=K’' with M= (4,17) 
A 4-t 

















where the rows of the 4-¢ by 4:3¢ matrix N consist of the norma- 
lized vectors n;. This transforms the 6-functions into (now explici- 
tely written) 6(K,’) .... 6(Kj..,). We further note that K? = Kj, 
as M is only non-orthogonal in S. Thus the integral (4,16) is trans- 
formed into (disregarding i’s) 
§(4-2t)(K’) 
“(3 tm?— K’2)3t ~ 
1 

(3 tm?— K%)3t 


(38)3¢|M|-1 [a%-29 K’ 


= |M|-*(80)* [aK 


f= Opt. 

= |M|-1(8#)34(4 2) gnnr (3tm?)-*. (4,18) 

For estimating |M| we may use Hadamards inequality which yields 

|M|? < 3*°?t, as there are 4-2t a’s with a? < 3 and n? = 1. This 
gives us the estimate 


K> (4 x)-*-# (8 1)" Co mt (4,19) 
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Comparing this expression with (4,3) we obtain 
# t \2¢ 
| ax |Q(a)|?> (gz) - (4,20) 


The first part of (4,14) can be estimated by observing that in view 
of the inequalities (4,9) the minimum value and the maximum 
value of Q with respect to x is given by 


min (2) a7 ( = 4) oe =") a= e) Q _ (1—zs) (%—2s) re (2— 23) ., (4,21) 


1 — 25) (2,—2_) — (2,— 24)?’ “max 


AS &pin 18 positive definite, the integral 


5 mt— py min —__ C(p%) m-* (4,92) 


{{Q’ |? (m? (1—a) + uy Uy (m?— p?)"/2+2} max 





[ dy dz du {u, a? (m?— p*)} min 
is positive for m? > p®. The integrand is independent of y and, 
therefore, [ dy can be replaced by 11). The remaining integral does 


not contain n nor is it dependent on the structure of the graph. This 
casts the estimate for H into its final form 


H> (. y" C (p2) et -2 3-Sn/2+11 (=)= 


\n 


ye (p2) A" n-2g 


Now we know that there are n/2—2 factorial graphs with the 
lower bound H. Therefore the series A, diverges at least”) as 


n2 


'S " cy (n= 2) . 


4amA 
The series does not converge for any value of 4. If any solution 
exists at all it must be nonanalytic for A= 0. On the other hand 
we know that for 4=0 there is a solution for the equations 
namely the free field. If there is any solution for the 4+0 it seems, 
therefore, doubtful whether it has any relation to the solution for 
4=0. It is to be expected that if there are solutions for 2+ 0 
that they have quite an unphysical behaviour. 


1) One can estimate easily that, for instance, C(o) > 1/5040 and C(p*) > C(o) 
for p? > 0. 

2) Hurst has found an upper bound for the terms of the series which shows 
that the series diverges actually in this manner. 
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It remains to be discussed whether the series converges in the 
domain p? > m*. One might suspect that if the series converges 
for p? > m? one could extend this expression to the region p? < m? 
by analytic continuation with respect to p. Though the author has 
not been able to prove anything in the case considered here, one 
can construct examples for which the series diverges for all values 
of the external momenta. If one investigates a theory with the non- 
linearity Ay* the convergence of the k-integration of an expression 
is only determined by the number of external momenta and not by 
the order of the graph. The present analysis can be carried through 
in this theory along similar lines, although things are more com- 
plicated in this case. But perturbation theory diverges for this 
theory on account of the same reasons. If one considers a graph 
with eight external lines, as for example 





ky bee 











h\ bi thyrky 


Fig. 4. 


he gets an expression of the form J da (md; ;(x) k,k;)-®. This is posi- 
tive for all values of the k’s and the same holds for contributions of 
higher order. The series diverges for the real part of this matrix ele- 
ment (which corresponds to the principal value of the above integral) 
irrespective of the value of k. But the real and the imaginary part 
of matrix elements are connected by certain equations and the series 
will diverge for the imaginary part too. 


To sum up one can say that the chances for quantized fields to 
become a mathematical consistent theory are rather slender. 


This work was carried out partly under a Fellowship of UNESCO 
and partly under a fellowship of the Swiss Federal Institute of 
Technology. The author wants to express his thanks to Prof. Pauni 
for his hospitality in his Institute and for making the work of 
Horst available to him. The author is very much indepted to Dr. 
KA.uEN for pointing out a mistake in the original manuscript. 
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Appendix a. 


In this appendix we shall show that the Eigenvalues £; and « 
(3,3) are positive. Considering p as one of the k’s we prove this by 
showing that the quadratic form of the k’s in the denominator is 
positive definite. The quadratic form is obtained by combining the 
denominators 


Ag (1 — ay) + dy(%1— 2p) + +++ + dy ty =k, kj Q5; (a,1) 


where d; = (Xk; c,)? and the c are 0,1. Q is a real symmetrical 
quadratic form and has, therefore, real Eigenvectors. If «; is the 
normalized Eigenvector belonging to the Eigenvalue #, 6B equals 
x,x;Q;;. If one takes for k the four k,; vectors (% , all having time 


direction we can write B =k;k;Q,;. For these k’s all the d’s in (a,1) 
are positive and f is bound to be positive as all the terms on the 
left hand side of (a,1) are positive. 

For proving 0<«(a#) <1 it suffices to observe that 0<s(#) <1 
and r;Qi.1; > 0 as we have shown above that a(x) > 0. The first 
of the above statements follows from the fact that s is of the form 
Xi (%q,—%a,11) With a; > a, for i > k. Keeping in mind that Q is 
positive definite we have in the diagonal system of Q:7;Q5) 1. = 
r;2 B>1 > 0 which is our second assertion. 


Appendix b. 


Here we shall investigate the case where in the k-integration some 
k remain in the numerator. Differenciating the general formula for 
k-integration with respect to a one gets 


— akkt —2in® 
[L-ak?}" ~~ a Lt-8(r—1) (r—2) (r—3) * 





(b, 1) 
Now we assume for simplicity that only one k occurs in the numer- 
ator and for this k we have a propagation factor of the form 


__ make 
[B m?— «a ke} * 


If one combines this line with the other lines of the graph we arrive 
at an expression of the form, 


d k(m?— k2) 
/ [bm?—(kQk+2pkr+e ple (b, 2) 
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We use matrix notation and write kQk for k;Q,;k;. By the sub- 
stitution k, > k,—Q7! rp we get rid of pk; in the denominator 
and a second substitution for all the other k’s and an unitary 
transformation U makes the denominator diagonal. The numerator 
becomes in this process 
2__}2 [72 " \? 12 

m*—k,Un—\@_) PB» 
if we drop linear in the k as they vanish by integrating. Applying 
the formula (b, 1) we are left with 


1 1 1 2 ry 2 2 bm?—«p* (U,,)? - 
Im —(a)P a ca oe (b, 3) 


The # are the Eigenvalues of Q. Now we know the inequality Q,;(2) 
> 1,(x). It can be seen immediately if one introduces the Feynman 
variables in the way similar (8, 1). Then Q,, is the total length of 
all lines containing k, while r, is the length of the lines containing r 
and k,. As the £ are positive the square bracket in (b, 3) is positive 
for p? < m?, 

The reason for the signs in favour of our general condition lies 
in the fact that we have a sign in (b, 1) opposite to the sign of an 
integral where no k? occurs in the numerator. 


If there are two integration variables in the numerator, (m?—k?) 
(m?—k?2) the term of the fourth order too yields the correct sign. 
This can be seen by the following analysis. The real unitary trans- 
formation transforms the term of fourth order in k into a term 
(k,k,,) (k,k,,) Uy ., Uy, Ug, Uz,-.. As the denominator then depends 
on the k squared only we make the substitution 


a 7.% TB ] a2 Te 1 1 a ae Py 
ks ki, ke kf —> Oss: Ops tae 8 ? 86 Oy 9s kt ke + 


+7 FP 28 9,5 Digs WER. 


Then the numerator of this term becomes 


aa0 ° , 1 
ke ks Ui, UE, + 2 U;, Ui, U;, U,,| " 


The square bracket is again positive as it may be written in the 
form 


2 ; 1 . 
E Ui, U+> (U,,U,,+ U,, U,,)? ~ 
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In the same way he can convince himself that in more compli- 
cated cases the signs of the terms are so as to make the expression 
positive for p? < m?. 

A more general argument runs as follows. We replace the propa- 
gation factor between self-energy graphs by 


1 
(m?— Ky ° 
‘ 
1°) 3 fe) 
and consider the expression as a function of the r. For r,2 <r < oo, 
there are no k? in the numerator and the expression can be calcu- 
lated in the same way as in (3,2). 


The expression thus obtained is positive for all values of r and 
goes to infinity for those values of r for which the k-integration 
becomes divergent. As this formal expression is an analytic function 
in the r which is identical with our expression for all integral r bet- 
ween 2 and oo, it has to be identical with our expression for r = 1. 





Construction de compteurs Maze et leurs caractéristiques 
de 0° 4 50° C. 
par R. Favre et Ch. Haenny (Lausanne). 
(22. XI. 1952.) 


Sommaire. Les tubes compteurs de verre 4 cathode externe dont l’anode est 
liée galvaniquement a son amarre de maniére a permettre un dégazage par la 
chaleur, donnent des décharges parasites qui ont été supprimées par la métallisation 
interne de l’extrémité correspondante de la coque. Ces compteurs sont caractérisés 
par des seuils identiques et des paliers de 300 volts environ dont la pente moyenne 


est respectivement de 3% et de 1,6% pour des résistances anodiques de 1 M et 
20 M. Nous donnons une explication des faibles variations de la pente, observées 
lorsque change le taux de comptage ou la température. En 10° impulsions, le 
palier se raccourcit au point de disparaitre. 


Introduction. 


Depuis plusieurs années*) que nous construisons en série des 
tubes compteurs en verre 4 cathode externe du type Mazex, nous 
avons éprouvé le besoin de les mieux connaitre et d’améliorer la 
qualité de ceux qu’au début, nous réalisions sur la seule base des 
données publiées. Il nous fallait pour des études sur la radiation 
cosmique, des compteurs présentant des seuils identiques et des 
paliers de faibles pentes dont les caractéristiques soient relative- 
ment stables dans le temps et dans un domaine de température 
compris entre 10° C et 40° C. 

Nous donnons ci-dessous les principes et les détails de construc- 
tion qui nous ont permis d’atteindre ce but, ainsi que les valeurs 
numériques intéressant ces compteurs. 


I. Effets de bouts et géométrie électrostatique. 


A cété de la variation du volume sensible des extrémités des 
compteurs, nous avons mis en évidence un effet de bout, consistant 
en impulsions parasites générées & proximité de l’amarre dans les 
circonstances précisées ci-dessous. 

La jonction métallique entre l’anode du compteur et |l’amarre 
(fig. 1) est justifiée par le fait qu’une matiére isolante organique 


*) Les premiéres recherches ont été effectuées en collaboration avec M. DEGAL- 
LIER, Helv. Phys. Acta. 22, 621 (1949). 
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(nylon) ne résisterait pas & la température élevée nécessaire pour 
un bon dégazage. Elle présente cependant un inconvénient, en ce 
que la résistance de contact entre l’amarre et l’enveloppe de verre 
est en général grande par rapport 4 la résistance entre ce point de 
contact et la cathode du compteur. La différence de potentiel appli- 
quée au compteur se concentre donc essentiellement entre l’amarre 
et la partie avoisinante de l’enveloppe. Le champ électrique intense 


Anode Ni Copper. Amarre Cu 











Fig. 1. 
Amarrage de l’anode. 


qui résulte de cette géométrie électrostatique particuliére est a 
Vorigine de décharges parasites dont le nombre augmente avec la 
tension. 

Ce phénoméne peut étre évité en appliquant la méthode proposée 
par D. Buanc et M. Scutrer en vue de mieux définir le volume 
sensible!) consistant & porter les extrémités extérieures des comp- 


teurs au potentiel du fil. Le méme résultat est obtenu plus simple- 
ment en réduisant la résistance de contact par ramollissement, sous 
vide, du verre en contact avec l’amarre ou par une métallisation de 
la partie interne du compteur en contact avec celle-ci (fig. 2). 


Région métallisée \ 











Fig. 2. 
Région métallisée. 


Le ramollissement du verre au voisinage de |’amarre donne d’ex- 
cellents résultats, mais entraine fréquemment un relachement du fil 
qui doit étre tendu 4 nouveau. 

Parmi les divers procédés de métallisation expérimentés, nous 
avons retenu l’argentage chimique, assurant le maximum de simpli- 
cité et de qualité du dépét. L’amélioration de fonctionnement des 
compteurs métallisés compense largement la faible augmentation 
du prix de revient résultant de cette opération supplémentaire dont 
le détail est donné par la suite. 
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II. Données constructives et mode opératoire. 
1, Préparation des anodes. 


Les anodes sont constituées d’un fil de molybdéne de 0,1 mm de 
diamétre, relié aux deux extrémités 4 un fil de copperclad de 0,3 mm 
de diamétre, dont |’un se termine par l’amarre et |’autre sort du 
compteur pour assurer la liaison électrique extérieure. La jonction 
entre le fil de molybdéne et ceux de copperclad est électrostatique- 
ment blindée par un cylindre de nickel écrasé d’un coup de pince 
dans sa partie médiane. 

L’extrémité des cylindres de nickel tournée vers |’intérieur du 
compteur est préalablement repoussée sur un tour rapide de ma- 
niére & constituer une calotte sphérique, percée d’un trou juste suf- 
fisant au passage du fil. 

Les anodes sont dégraissées par immersion dans le benzéne ou 
le tétrachlorure de carbone qui doit étre complétement éliminé 
avant d’introduire les anodes dans les enveloppes de compteurs. 


2. Enveloppes des compteurs. 


Les enveloppes doivent étre de faible résistivité électrique. Le 
verre « Novo» des Cristalleries et Verreries de Choisy le Roy (France) 
déja utilisé par Mazz est satisfaisant & ce point de vue et en raison 
de sa parfaite adhésion au copperclad. Le soufflage des enveloppes 





a 





Fig. 3. 
Coque de compteur avant son montage. 


est précédé d’un lavage du verre a l’eau courante pour éliminer les 
poussiéres. Les coques, qui présentent la forme indiquée par la 
fig. 8, sont ouvertes aux deux extrémités, l’une d’elles se terminant 
par un tube de quelques centimétres. 


3. Nettoyage des enveloppes. 


Ainsi que l’a indiqué Mazx, l’importance du nettoyage des en- 
veloppes de verre neuves ne doit pas étre exagérée. Nous avons 
construit d’excellentes séries de compteurs, en ayant rincé les en- 





56 R. Favre et Ch. Haenny. 


veloppes une seule fois 4 l’eau distillée. I] est prudent cependant 
d’utiliser un agent de nettoyage plus énergique, tel que le mélange 
chromique, dans lequel les compteurs séjourneront quelques heures. 
Cette mesure s’avére indispensable s’il s’agit d’enveloppes provenant 
de compteurs usagés ou encore en vue d’un argentage ultérieur. 


Le nettoyage 4 l’acide chromique doit étre suivi d’un ringage 
abondant a l’eau courante, puis 4 l’eau distillée. Les enveloppes 
sont alors séchées, si possible dans une étuve, a l’abri des pous- 
siéres. 


4, Argentage de la zone de contact avec l'amarre. 


L’argentage est précédé d’un paraffinage du tube terminal de 
l’enveloppe, opération destinée a délimiter l’une des frontiéres de 
la partie argentée. La solution d’argentage est alors aspirée dans 
les coques fixées verticalement, jusqu’a la hauteur désirée. La durée 
du dépét est d’environ 10 minutes aprés lesquelles la solution est 
éliminée et les enveloppes rincées & l’eau bouillante qui entraine 
la paraffine, puis 4 l’eau distillée. Les enveloppes argentées sont 
& nouveau séchées avant de recevoir les anodes. 


La solution d’argentage?) est constituée par le mélange, au der- 
nier moment, de volumes égaux des deux solutions suivantes: 


a) 10 cm? d’une solution d’AgNO, 4 10%, 


30 cm? d’eau distillée, 

NH,OH jusqu’a redissolution presque compléte du précipité, 
25 cm d’une solution de KOH a 10%, 

NH,0H jusqu’a redissolution presque compléte du précipité. 
Filtrer la solution. 


Solution réductrice comprenant: 


10 cm* d’une solution de glucose 4 5%, 
12 cm? d’une solution de sucre de canne 4 10%, 
43 cm? d’eau distillée. 


5. Montage de l’anode. 


Les anodes sont introduites dans les enveloppes et scellées tout 
en exergant une traction qui les tend suivant l’axe de la coque. 
Une gorge destinée a faciliter le scellement des compteurs aprés rem- 
plissage est pratiquée a l’endroit convenable du tube terminal. 
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6. Dégazage. 


y 


Les enveloppes sont soudées & la rampe 4 vide par le tube ter- 
minal resté ouvert, vidées & une pression de 10-4 mm de Hg, puis 
chauffées pour assurer le dégazage des parois. A défaut d’un four 
spécialement destiné a cet usage le dégazage peut étre effectué par 
le chauffage au chalumeau de chacun des compteurs pendant 10 mi- 
nutes, jusqu’a légére coloration orange de la flamme initialement 
bleue. Le pompage doit étre maintenu durant le dégazage. 

Un compteur mal dégazé est rapidement mis’ hors d’usage par 
un taux de comptage élevé & une température supérieure a 30° C. 


7. Remplissage. 


Le remplissage par introduction successive dans les compteurs 
de Valcool éthylique et de l’argon est a déconseiller, la proportion 
des gaz variant de l’un a l’autre selon leur position sur la rampe. 
Il est préférable de préparer le mélange des gaz en proportion con- 
venable (1 volume d’alcool éthylique pour 9 volumes d’argon) dans 
un récipient suffisamment volumineux pour permettre le remplis- 
sage des compteurs 4 la pression de 10 cm de Hg environ. 

Nous avons utilisé avec succés l’argon technique pour soudure, 
dune pureté de 99,8% et l’alcool éthylique absolu, vendu dans le 
commerce en ampoules scellées. 


8. Derniéres opérations. 


Les compteurs remplis peuvent étre scellés et du méme coup dé- 
tachés de la rampe. Leur surface extérieure est soigneusement net- 
toyée et les extrémités déséchées au chalumeau sont plongées dans 
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Copperclad Ozokerite Graphite Molybdéne Czokérite Argentage 
Fig. 4. 
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Type de compteur Maze & cathode externe. 





l’ozokérite en fusion. La cathode est constituée par l’application 
d’une couche conductrice de graphite sous forme «d’aleool DAG». 
Les compteurs préparés selon ce procédé (fig. 4) sont utilisables 
sans délai. 
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Ill. Caractéristiques des compteurs. 


1. Courbes de comptage et durée de vie. 


Les mesures de palier données ci-dessous ont été faites sur 5 comp- 
teurs non sélectionnés, avec une échelle de 64 précédée d’un préam- 
plificateur. La résistance anodique R (fig. 5) a été successivement 
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GM : tl 
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AT 
Fig. 5. 
choisie de 1 M et 20 M puis remplacée par un circuit coupeur élec- 
tronique (fig. 6). 
Les résultats obtenus 4 1’effet zéro et & un taux de comptage d’en- 
viron 7000 imp/min. sont résumés dans le tableau 1. 
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Deux causes peuvent expliquer la réduction de pente qui, pour 
une résistance anodique de 1 M, caractérise |’élévation du taux de 
comptage. 

La premiére réside dans ]’existence d’impulsions parasites dont 
la proportion diminue lorsque le taux de comptage augmente. La 
présence de ces impulsions parasites est confirmée par |’améliora- 
tion de la pente, qu’entraine l’usage d’une résistance anodique 
élevée ou d’un dispositif coupeur électronique. 











Fig. 6. 
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La seconde cause, comme nous le verrons dans le paragraphe 
suivant, est & rechercher dans la présence de la coque de verre sur 
le chemin des ions. Les compteurs en sont affectés aussi bien pour 
une résistance anodique de 1 M que de 20 M. §8’il ne semble pas en 
étre ainsi dans ce dernier cas, c’est que le temps d’inefficacité da a 
Vimpulsion négative du potentiel d’anode est alors plus grand que 
le temps mort du compteur et diminue avec la surtension, ce qui 
a pour conséquence d’accroitre la pente lors d’un taux de comp- 
tage élevé. 

Tableau 1. 





Longueur 


Circuit | Imp/min. | ape Pente du palier 


i 


R= 1M 300 20° C 3,1% + 0,2 250 V 
R= 1M 7000 20° C 2,2% + 0,1 250 V 
R=20M 300 20° C 1,6% + 0,2 300 V 
R=20M 20° C 1,8% 300 V 
R=20M 40°C 12 % 300 V 
R=20M 300 55° C 1,5% 300 V 

coupeur 300 20° C 22% 150 V 

coupeur | 7000 20°C | 2 MS 150 V 











La longueur du palier obtenu avec le dispositif coupeur électro- 
nique est liée 4 l’amplitude de l’impulsion négative de coupure. 
La vie des compteurs dépend dans une large mesure de la sur- 
tension de fonctionnement et du circuit d’alimentation (valeur de 


V/A Tension dy seuil 





510° 10° ump. recues 
Fig. 7. 


Elévation du seuil avec l’usure du compteur. 


la résistance anodique ou dispositif coupeur électrique). L’usure se 
manifeste par une élévation du seuil, la partie supérieure du palier 
n’étant que peu affectée. Il en résulte une réduction de la longueur 
de celui-ci. La figure 7 donne la courbe d’élévation du seuil en fonc- 
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tion du nombre d’impulsions regues, pour des compteurs travaillant 
avec une résistance anodique de 1 M et une surtension maintenue 
a 100 V. Dans ces conditions le palier proprement dit disparait aprés 
10° impulsions environ. Le renouvellement du gaz d’un compteur 
usé n’entraine pas sa régénération. I] est nécessaire pour cela de 
nettoyer & nouveau l’enveloppe a l’acide chromique. La méme 
constatation a été faite par D. Bhanc pour des compteurs d’un 
type semblable, contenant du méthane au lieu d’alcool. 


2. Influence du taux de comptage sur la caractéristique. 


Ainsi que R. Mazx) l’indique dans le premier mémoire sur ce 
sujet, un compteur a parois de verre et cathode externe se diffé- 
rencie des autres compteurs habituels par la présence, du cété de 
la cathode, d’une résistance élevée due a l’enveloppe de verre. 
La paroi interne de la coque se charge lors du fonctionnement du 
compteur et il y a lieu de faire intervenir en plus de la résistance 
une capacité en paralléle avec celle-ci. Nous avons mesuré ces gran- 
deurs comme |’a fait Maze en remplissant le tube d’une solution 
électrolytique conductrice de |’électricité. Nos valeurs s’accordent 
avec celles de cet auteur. Pour un compteur de 45 cm de longueur 
et de 2,8 cm de diamétre: 


Capacité C = 3500—4000 pF 


Résistance transversale R = 50 mégohms a 20° C 


ce qui correspond a une constante de temps de 0,2 sec. dont il peut 
étre nécessaire de tenir compte dans l’utilisation et l’interprétation 
des résultats expérimentaux. 


Pour un taux de comptage n de plusieurs dizaines d’impulsions : 
sec., on peut assimiler l’écoulement des charges au travers de l’en- 
veloppe de verre & un courant continu 1, d’ot la relation suivante 
entre la tension V, appliquée au compteur et la tension effective V, 
existant entre la paroi interne de la coque et le fil axial: 


V.=V,— Ri. (I) 


Nous avons mesuré le courant de décharge en placant un galvano- 
métre entre la cathode et la terre. La modulation du courant, due 
aux impulsions successives, est soigneusement filtrée. La mesure pré- 
cise du courant peut se faire pour un taux de comptage supérieur 
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& 1500 impulsions/min. Les résultats obtenus sont indiqués par la 
fig. 8. On voit que le courant de décharge a taux constant est appro- 


72000 imp /min 


Pe 


4000 imp mn 


1600 /mp/min 








1100 100 ib 


Fig. 8. 
Courant de décharge d’un compteur de 30 cm. 


ximativement proportionnel 4 la surtension, sur la plus grande 
partie du palier. 
i=u(Vi—V,). (II) 
Vy = seuil de la région de GEIGER; 
u = coefficient de proportionnalité pour la premiére partie de la courbe. 


On retrouve le point anguleux observé par d’autres, FENTON et 
Fuuier‘). La tension de ce point anguleux s’éléve avec le taux de 
comptage. 

de (I) et (ID) 
Ve = Va (1 — uR) + uRVy 


¥.=AV,+B 
La fonction n(V,) correspond aux indications que fournirait un 


compteur a cathode interne, toutes choses restant égales par ailleurs. 
Le paramétre «A» tend vers 1 et B vers 0 pour un taux de comptage 
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tres faible (effet zéro) ou une faible résistance de la coque de verre 
(température élevée). Dans ces deux cas, le compteur se comporte 
du point de vue électrique comme un compteur & cathode interne. 
Dans les autres cas, «A» est inférieur 4 l’unité. I] résulte de la re- 
lation d’affinité qui lie V, et V, que la longueur du palier augmente 
avec le taux de comptage tandis que la pente diminue. 

Il est bien des cas pour lesquels cette déformation est sans im- 
portance et ne nécessite aucune correction. I] n’en est plus de méme 
sila mesure doit s’effectuer & plusieurs taux de comptage trés diffé- 
rents, comme cela se présente lors de la détermination d’une période 
d’un élément radioactif & décroissance rapide. La valeur de V, varie 
du début a la fin de la mesure. 

Les résultats ne seront valables que si le palier est assez long 
et de pente trés faible. Dans tous les autres cas il est nécessaire de 
corriger les valeurs trouvées a l’aide des relations ci-dessus pour 
ramener chaque résultat 4 celui que |’on aurait obtenu pour une 
tension V, constante. Cette correction bien que facile est fastidieuse 
et ne peut qu’agrandir la marge d’erreur. 
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Fig. 9. 
Résistance transversale de la coque. 


3. Effet de température. 


L’influence de la température peut se décomposer en un effet 
systématique dti 4 la variation de résistance de la coque (fig. 9) 
qui ne se manifeste que pour un taux de comptage élevé et en un 
effet comprenant les autres causes possibles. 
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Le comportement des compteurs aux basses températures a été 
étudié en les plagant dans une enceinte frigorifique dont la tempé- 
rature est réglable & + 1° C. 


Les paliers mesurés pour l’effet zéro ne sont pas affectés jusqu’a 
la température minimum de 5° C, comme l’ont observé d’autres 
expérimentateurs. Au-dessous de cette température la détérioration 
du compteur n’est pas instantanée et résulte d’une photosensibili- 
sation du compteur, probablement due a une adsorption massive 
des molécules d’alcool. Les compteurs détériorés par |’abaissement 
de la température se restaurent spontanément a la température de 
20° C, mais il leur faut plusieurs jours pour retrouver leurs caracté- 
ristiques initiales, tandis qu’il suffit de quelques heures si les comp- 
teurs sont placés dans une étuve a 50° C. 


~~ 3 


Seul l’effet systématique di a la variation de résistance électrique 
de l’enveloppe a été observé a chaud (tableau 1). 


Conclusion. 


Construits ainsi que nous |’indiquons, les compteurs & cathode 
externe du type Maze sont & méme de répondre aux exigences qui 
nous sont imposées dans notre étude du rayonnement cosmique. 
Le rendement de fabrication est voisin de 100% et les caractéris- 
tiques sont bien définies dans un domaine thermique de 5° a 50° C 
pour les compteurs 4 argon-alcool. Les précautions prises pour 
définir convenablement la géométrie électrostatique, permettent 
d’obtenir une pente assez faible pour ne pas justifier l’usage de gaz 
halogénes dont l’emploi entraine de sérieuses complications. Un 
dispositif coupeur électronique ou, si le temps mort n’est pas cri- 
tique, une résistance anodique élevée, peuvent réduire la pente 
et allonger le palier. Les faibles variations de pente résultant des 
changements de température ou du taux de comptage, s’expliquent 
toutes par la résistance électrique de la coque de verre et par l’effet 
d’impulsions parasites secondaires. 


La vie des compteurs dépend de la surtension et de leur géomé- 
trie. Les compteurs de 45 x 2,8 cm, travaillant a une surtension 
constante de 100 V, accusent en 3.108 impulsions, une réduction 
du palier de 100 V par élévation du seuil. C’est l’altération a laquelle 
on peut s’attendre lors de mesures sur la radiation cosmique, effec- 
tuées en permanence pendant un an 4 l’altitude du Jungfraujoch 
(3500 m). 
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Berechnung linearer, realisierbarer Netzwerke zur Erzielung 
optimaler Signal/Rauschverhaltnisse 
von K. Halbach, Physikalisches Institut der Universitit Basel. 
(20. I. 1953.) 


Zusammenfassung: Es wird eine allgemeine Methode zur Berechnung derjenigen 
linearen, realisierbaren Anordnung entwickelt, die am Ausgang ein optimales Ver- 
haltnis von Signal zu Untergrund ergibt. Die Eigenschaften von Signal und Unter- 
grund am Eingang werden dabei als bekannt angenommen, die Messzeit kann be- 
liebig vorgegeben werden. Als Anwendung wird der optimale Ionisationskammer- 
verstarker berechnet und mit dem RC—RC-Verstarker und dem delay line clipping 
verglichen. 


1. Problemstellung. 


Viele physikalische Messungen werden auf das Problem zuriick- 
gefiihrt, gewisse Eigenschaften eines elektrischen Signals, z. B. 
Kurvenform oder Signalhéhe, zu bestimmen. Da dem Signal immer 
ein Untergrund iiberlagert sein wird (z. B. Widerstandsrauschen, 
Réhrenrauschen), ist die gewiinschte Messung in den meisten Fallen 
nur mit begrenzter Genauigkeit méglich. Diese Messgenauigkeit 
wird im allgemeinen geandert, wenn man das vom Untergrund tiber- 
lagerte Signal durch ein Filter (Verstarker) gibt, und man kann 
nun diejenige lineare Anordnung*) suchen, welche die zu unter- 
suchende Signaleigenschaft mit optimaler Genauigkeit zu messen 
gestattet. Dabei kann man voraussetzen, dass die Anordnung selbst 
keinen zusaétzlichen Untergrund erzeugt. Die Berechnung der opti- 
malen Anordnung zur Bestimmung der Signalform wurde von Bopr 
und SHannon!) durchgefiihrt. Oft ist jedoch die Signalform be- 
kannt, und man wiinscht die Signalamplitude mit méglichst grosser 
Genauigkeit zu bestimmen. DEN Harroc und MU.uer?) haben 
dieses Problem bereits behandelt. Es soll hier erneut diskutiert 
werden, jedoch unter Beachtung zweier zusiétzlicher Bedingungen, 
welche die vom theoretischen und praktischen Gesichtspunkt aus 
unbefriedigenden Lésungen ausscheiden. 

Zunichst einmal muss man von der gesuchten Anordnung ver- 
langen, dass sie realisierbar ist; sie darf also kein Ausgangssignal 


*) Unter linearer Anordnung versteht man ein Netzwerk, das im Sinne der 
Fourierzerlegung durch eine komplexe Ubertragungsfunktion beschrieben werden 
kann, d.h. also, dass Eingangs- und Ausgangssignal miteinander durch lineare 
Gleichungen bzw. Differentialgleichungen verkniipft sind. 
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geben, das zeitlich vor dem erzeugenden Eingangssignal beginnt. 
Ferner muss eine praktisch brauchbare Schaltung natiirlich stabil 
sein, d. h. sie darf z. B. auf eine kleine Stérung nicht mit einem 
bis ins Unendliche andauernden, exponentiell anwachsenden Aus- 
gangssignal antworten. 

Weiterhin wird man in vielen Fallen mit grésserer Messzeit auch 
eine gréssere Messgenauigkeit erreichen kénnen. Man muss also, 
um nicht zu unendlich langen Messzeiten zu kommen, die Messzeit 
irgendwie beschranken. Das geschieht hier, indem die Zeit, zu der 
man die Signalhéhe misst, relativ zum zeitlichen Auftreten des 
Signals am Eingang beliebig, aber fest vorgegeben wird. 

Ist Sy (w) das Fourierspektrum des Signals I’, (t) am Eingang 
der gesuchten Anordnung, deren Ubertragungsfunktion () sei, 
so ist der zeitliche Verlauf des Signals am Ausgang des Verstiirkers 
gegeben durch: 


Hi(t)= | So(w) w(mje''dy (Qav=a). (1) 
Da der mittlere quadratische Fehler fiir die Messung der Signal- 
hdhe in einem vorgegebenen Zeitpunkt gleich dem Effektivwert des 
Untergrundes ist, wahlen wir als Mass fiir die Messgenauigkeit das 
Verhaltnis vom Quadrat der gemessenen Signalhéhe zum quadra- 
tischen Mittelwert (Quadrat des Effektivwertes) des Untergrundes. 
Beschreibt man das Rauschen am Eingang durch sein quadratisches 
Fourierspektrum N,(@), das wir als bekannt voraussetzen, so ist 
der quadratische Mittelwert des Rauschens am Ausgang bekannt- 
lich*) : 
R= | No(o) | u(w) |? (2) 


. 
—oo 


und man erhalt fiir das (quadratische) Verhaltnis von Signal zu 
Rauschen am Ausgang, wenn man in (1) fiir ¢ den zur Messung 
der Signalhéhe vorgegebenen Wert ¢, einsetzt: 


[ Sy() h(a) ci tody 2 

V H2(t ) so | 

*(to) = a _ '- | | 3 
| No(m)| w(o) [2 dv 


—c 


*) Aus formalen Griinden integrieren wir hier von — oo bis + oo und nicht, 
wie allgemein iiblich, von 0 bis + co. Dies bedeutet in der Normierung von N,(«) 
einen Faktor 2. 
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In diesem Ausdruck muss nun diejenige Funktion p(w) gefunden 
werden, die ihn zum Maximum macht. Dabei sind jedoch nicht alle 
denkbaren Funktionen zugelassen, sondern nur solche, die realisier- 
baren und stabilen Anordnungen entsprechen. 


2. Berechnung der Ubertragungseigenschaften des optimalen Netzwerkes. 


Die Behandlung des Problems wird stark vereinfacht und die 
Resultate werden sehr anschaulich, wenn wir die gesuchte Anord- 
nung durch zwei hintereinandergeschaltete Netzwerke ersetzen der- 
art, dass das quadratische Fourierspektrum des Rauschens nach 
der ersten Anordnung gleich eins, d.h. also frequenzunabhingig, 
wird (Fig. 1). Wir verlangen von diesen beiden Anordnungen, dass 
jede fiir sich realisierbar und stabil ist, so dass dies auch von der 
gesamten Anordnung gilt. Sind die Ubertragungsfunktionen der 
beiden Netzwerke y,(@) und ,(@), so ist: 


f4(@) = Hy (@) He (@) - (4) 


Das quadratische Fourierspektrum des Rauschens nach der ersten 
Anordnung ist gegeben durch: 


N,(@) = No(@) | 4, (@) |?. 


Da wir dies zu eins machen wollen, muss fiir “,(@) auf der reellen 
Frequenzachse gelten: 


| 4, (@) |? = 1/Ng (a) . (5) 


Betrachten wir die Ubertragungsfunktion y,() als Funktion der 
komplexen Variablen w, so ist die Forderung der Realisierbarkeit 
und Stabilitét bekanntlich aquivalent der Bedingung, dass “,(@) 
in der unteren w-Halbebene (d. h. fiir w mit negativem Imaginirteil) 
keine Pole hat. Ebenso miissen wir aber auch von ,(@) verlangen, 
dass es keine Nullstellen in der unteren w-Halbebene hat; diese 
Nullstellen kénnten naémlich durch y.(@) nicht mehr kompensiert 
werden (die Anordnung mit der Ubertragungsfunktion s.(@) wiire 
dann instabil oder nicht realisierbar) und wir wissen nicht, ob die 
optimale Ubertragungsfunktion m(w) = 44(@) M(m) Nullstellen in 
der unteren o-Halbebene haben darf. Durch Gleichung (5) und 
die Bedingung, dass ,(@) in der unteren w-Halbebene keine Null- 
stellen und Pole haben darf, ist ,(@) bestimmt. In den meisten 
Fallen ist No(w) der Quotient zweier Polynome von w?, so dass 
man dann /4,(@) sehr leicht direkt angeben kann. In komplizierteren 
Fallen lisst sich y,(@) mit Hilfe der Beziehungen zwischen Ampli- 
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tudengang und Phasengang, die man z. B. bei Bove’) findet, 
berechnen. 

Durch das Netzwerk y,(m) wird das Eingangssignal F’)(t) verformt; 
am Punkt 1 (Fig. 1) gilt namlich fiir das Fourierspektrum des 
Signals S,(w) = So(m)-,(@), und fiir das Signal selbst ist: 


F(t) = [So(o) (0) ay. (6) 


f(a) ie poy 


' lo 
a Vitg)=fhtede 
Heo 0 eo ie 
; 5, (0) =S0 ht pel 4 f1G,/2 
1 =f 2! 
Fi aes wate f oI ov 









































Fig. 1. 
Schematische Darstellung der allgemeinen Resultate. 


Den zweiten Teil des optimalen Verstirkers beschreiben wir zweck- 
miassigerweise nicht durch seine Ubertragungsfunktion .(@), son- 
dern durch deren Fouriertransformierte: 


oo 


f(t) = | pa(w)e'e'dr. (7a) 
Diese Funktion hat eine direkt anschauliche Bedeutung: da ein 
zur Zeit Null erfolgender 6-Stoss (unendlich kurz dauernder Impuls 
mit endlicher Flache) ein frequenzunabhingiges Fourierspektrum 
hat, ist f(t) das Ausgangssignal, wenn man auf den Eingang der 
Anordnung mit der Ubertragungsfunktion p(w) zur Zeit Null einen 
6-Stoss mit der Flache 1 gibt. Aus dieser Deutung wird sofort klar, 
was die Realisierbarkeit fiir f(t) bedeutet: fiir t << 0 miissen wir 
f(t) = 0 setzen. Die Stabilitatsbedingung ist ebenso leicht einzu- 
sehen: fiir t+ co muss f(t) gegen einen endlichen Wert gehen; 
wie wir spiter am Resultat sehen werden, wird diese Bedingung 
immer von selbst erfiillt sein. Wegen der Realisierbarkeit erhalten 
wir fiir die Riicktransformation von (7a): 


Hao) — [Feri dt. 


Mit Hilfe dieser Beziehungen werden wir nun den Ausdruck fiir 
H(t) (Gleichung (1)) umformen: 


H (t) = / So (w) p(w) ee'dy= [Soi mge*tdy. 
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Hierin ersetzen wir Ms durch (7b): 


H(t) = [ [8(0) ow) ty (@) f (x) eb dadr. 
«& 0 
Vertauschung der Integrationsfolge und Durchfiihrung der Inte- 
gration iiber » ergibt mit (6): 
oo t 
H(t) = [F,t—2)f(2)de=/F,(a)f(t—a)de. —®) 
0 —oo 
Auch diese Gleichung erlaubt wieder eine direkte Interpretation: 
denkt man sich F’,(x) aus 6-Stéssen der Flache F(x) -dx zusammen- 
gesetzt, so gibt ein derartiger, am Eingang zur Zeit « erfolgender 
Impuls am Ausgang zur Zeit ¢ den Beitrag F',(x)-f(t — x)-dz, und 
Gleichung (8) ist die Summe aller Beitrage der gesamten Vergangen- 
heit. 
Um den Ausdruck fiir das Rauschen (Gleichung (2)) umzuformen, 
verfahren wir ganz ahnlich: mit (2) und (4) ist: 


B= | No(o) | 144) |? | s29(e) | dr. 
Wegen (5) wird: . 
R= | |pa(o) 2dv= | palo) ta(—0) dy. 


ford 


Wenden wir hierin auf (— w) Gleichung (7b) an, so wird: 
_ 


R=| [f(a 2) fy (w) e'”* da dy. 


—o U 


Vertauschen wir darin die Integrationsfolge und beachten bei der 
Durchfiihrung der Integration iiber » Gleichung (7a), so erhalten 
wir: > 
R= / f2(a) da. (10) 
0 
Setzen wir (8) und (10) in (3) ein, so bekommen wir fiir das (quadra- 
tische) Verhiltnis von Signal zu Rauschen: 


; [Frau-arerae| 
acta eal : 
| P(e) de 
0 
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Um diesen Ausdruck durch geeignete Wahl von (2) zum Maximum 
zu machen, variieren wir in der iiblichen Weise: wir ersetzen darin 
f(x) durch f(x) + eg(x), differenzieren V? nach e, lassen ¢ gegen 
Null gehen und setzen dann das Ganze gleich Null. Wir erhalten so: 
, 1H , 

(H.-5 %-°B)_-9- (11) 

Mit (8) und (10) wird: 
(H’).-0= | Fi(ty—a) g(#) da ; (R,),-9=2 | f(a) g() da. 
0 a > 

Setzen wir diese beiden Ausdriicke in (11) ein und schreiben alles 
unter ein Integralzeichen, so muss sein: 


/ [Fs (o—2)— nt (2)| g(a) da=0. (12) 


Hierin ist H/R beziiglich der Integrationsvariablen « konstant und 
enthalt die Funktion g(x) nicht mehr. Eine in f(t) noch enthaltene 
. willkiirliche Normierungskonstante kénnen wir daher so festlegen, 


dass H/R=1 (13) 


wird. Da Gleichung (12) fiir jede Funktion g(a) erfiillt sein muss, 
bekommen wir als gesuchte Funktion, wenn wir wieder ¢ statt x 


schreiben: {(t)=F,(t—t) (t>0). (14) 


Wegen der Realisierbarkeit ist, wie schon weiter vorne angegeben, 
f(t) =0 fiir t< 0; die Stabilititsbedingung ist auch erfiillt, da jede 
sinnvolle Funktion I’,(t) fiir sehr grosse negative Werte von ¢ ver- 
schwinden wird. Aus (14) lasst sich mit Hilfe der Beziehung (7b) 
die Ubertragungsfunktion p(w) ohne Schwierigkeiten berechnen. 
Wir verzichten hier auf Herleitung und Angabe des Ausdruckes 
fiir .(@), da er nichts wesentlich Neues erkennen lasst und, im 
Gegensatz zu der sehr anschaulichen Gleichung (14), sehr uniiber- 
sichtlich ist. 

Nachdem jetzt die Ubertragungseigenschaften des gesuchten 
Systems bekannt sind, miissen wir noch bestimmen, wie das 
damit erzielte Verhiltnis von Signal zu Rauschen abhangt von 
der Zeit t), zu der man die Signalhéhe misst. Mit den Glei- 
chungen (8), (8), (13) und (14) erhaélt man dafiir: 


V (to) = |/ ' (15) 
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Es sei an dieser Stelle noch einmal daran erinnert, dass, wie schon 
angedeutet, 1/V(t)) der mittlere relative Fehler fiir eine Messung 
der Signalhéhe ist. Wenn man den Grenzfall t) > co betrachtet, 
wird: 


V2(c0) 


Eine ganz gleiche Betrachtung, wie wir sie zur Herleitung 
Gleichung (10) aus Gleichung (9) angewandt haben, fiihrt zu: 


V2 (co) = / |S, (@) |? dv = / (So (@) |? [4a (@) |? dv. 

Wegen (5) wird somit: 
me / ris (w) |? . 
V(co) = | / Y, (0) dy. (16) 
Obgleich bei vielen praktischen Problemen die Gleichung (15) 
von besonderem Interesse sein wird, ist doch auch der Ausdruck 
(16) sehr bemerkenswert: er erméglicht es auf einfache Weise, 
direkt aus den vorgegebenen Gréssen So(w) und No(@) zu berechnen, 


wie gross das mit einer linearen Anordnung prinzipiell iiberhaupt 
erreichbare Verhialtnis von Signal zu Rauschen ist. Fig. 1 gibt noch 
einmal eine Ubersicht, wie sich die Gréssen p4(), f(t), V2(t)) und 
V2 (co) aus den vorgegebenen Funktionen F,(¢) und N,(w) berech- 
nen lassen. 


3. Anwendung auf einen Spezialfall. 


Zur Erliuterung der allgemeinen Methode soll im folgenden als 
einfaches Beispiel die Messung von Ladungen mit Hilfe der Ionisa- 
tionskammer behandelt werden. Im allgemeinen verbindet man 
bei derartigen Messungen die Sammelelektrode der Ionisations- 
kammer mit dem sonst freien Gitter der ersten Verstarkerréhre, 
wahrend die andere Elektrode an eine Gleichspannung gelegt wird 
und somit wechselstrommissig geerdet ist. Um eine etwas ein- 
fachere Darstellung zu bekommen, nehmen wir an, dass man von 
der Anode der ersten Réhre zunichst einmal auf einen idealen (d. h. 
frequenzunabhingigen) Breitbandverstirker BBV geht, so dass man 
an dessen Ausgang geniigend tiber dem Stérpegel der folgenden 
Anordnung liegt und von der ersten Réhre entkoppelt ist (Fig. 2). 

Wie Miuatz und Keuier‘) gezeigt haben, ist bis auf unwesent- 
liche Konstanten an der ersten Anode und damit auch am Ausgang 
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des BBV das quadratische Fourierspektrum des Rauschens gegeben 
durch: 
No(@) =1+a?/o?. (17) 


Unter der Annahme einer unendlich kurzen Sammelzeit der Kam- 
mer erzeugt ein zur Zeit t = 0 auftretendes ionisierendes Teilchen 
am Ausgang des BBV ein Signal von der Form: 


Fo(t) =1; (Fo(t) =0 fiir t<0). (18) 
Fiir den ersten Teil des optimalen Ubertragungssystems finden wir 
mit (5) und (17) sofort: . 
fy (@) =ia/(iw+a). (19) 






































Fig. 2. 


Ubersicht iiber den optimalen Ionisationskammerverstarker. 


Wie man sich leicht iiberlegt, entspricht das in Fig. 2 angegebene 
RC-Glied gerade dieser Ubertragungsfunktion. Das durch (18) 
gegebene Signal wird durch dieses RC-Glied bekanntlich umge- 
formt in: 

F,()=e" ; (F,() =0 fir t<0). (20) 
Damit und mit (14) ist nun auch das gesuchte f(¢) bekannt: 


f@=e*-e* (O<t<h) 
f() =0 (<0; >t). 
In Figur 3 ist diese Funktion fiir fj = 1/a dargestellt. Durch Ein- 
setzen von (20) in (15) erhalt man fiir das mit der optimalen Anord- 
nung erreichte Verhaltnis von Signal zu Rauschen: 
V2 (t,) = (1—e-2"") 2a. (21) 


Zum Vergleich mit in der Praxis verwendeten Verstarkern ziehen 
wir zunichst den sogenannten RC—RC-Verstirker heran. Bei ihm 
wird der Frequenzgang nach hohen und tiefen Frequenzen je mit 
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einem RC-Glied beschnitten; seine Ubertragungsfunktion hat also 
die Gestalt: 


1(@) po = q+ t wf (i wo + 4) (VO + we). 


Die beiden Kreisfrequenzen ,, w, sind die reziproken Werte der 
beiden verwendeten Zeitkonstanten R, C, baw. R,C,. Wie van 
HEERDEN®) gezeigt hat, wird das Verhiltnis von Signal zu Rauschen 
am gréssten, wenn man @, = w, = 4 macht. In diesem Fall hat 
das Ausgangssignal seine Maximalamplitude zur Zeit t = 1/a, und 
es wird: 
Vig = 2/a6*. 

Durch Vergleich mit (21) sieht man, dass fiir tf) > 1/a der optimale 
Verstaérker um den Faktor V,,./ Vee = e/2, d. hh. rund 36% besser 


fd) 
1+. 














a5 / 
Fig. 3. 
f(t) fiir den optimalen Verstarker (————) 
fiir den RC—RC-Verstarker 
fiir delay line clipping 
(t= 1 / a). 


ist als der RC—RC-Verstirker. In seinen Betrachtungen tiber die 
prinzipielle Grenze der Messgenauigkeit, wenn Signal und Rauschen 
durch (17) und (18) gegeben sind, kommt van HererpEN5) zu dem 
gleichen Resultat. Schreibt man fiir den optimalen Verstirker die 
gleiche Messzeit ty = 1/a vor, wie sie der beste RC—RC-Verstarker 
hat, so ist der Gewinn im Verhiltnis von Signal zu Untergrund nur 
noch 26%. Da das die tiefen Frequenzen beschneidende RC-Glied 
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des besten RC—RC-Verstiirkers identisch ist mit dem Gleichung (19) 
darstellenden RC-Glied des optimalen Verstirkers, kénnen wir den 
zweiten Teil der Ubertragungsfunktion des RC—RC-Verstirkers 
(= a/(a+i@)) mit dem f(t) des optimalen Verstirkers vergleichen. 
Die a/(a+iw) aquivalente Funktion f(t) ist in Fig.3 wiedergegeben. 
Man sieht aus dieser Darstellung, dass man als niichst bessere Appro- 
ximation des optimalen Verstarkers ein f(t) wahlen wird, das etwa 
der gestrichelten Geraden entspricht. 

Wie man leicht findet, ist eine derartige Anordnung, zusammen 
mit dem durch (19) beschriebenen RC-Glied, gleichbedeutend mit 
dem in der Praxis schon lange gebrauchten delay line clipping. Das 
damit im Maximum (bei ft) = 1,25/a) erreichbare Verhaltnis von 
Signal zu Untergrund liegt nur noch 10% unter demjenigen des 
optimalen Verstirkers bei unendlich langer Messzeit. Wegen der 
Verbesserung der Statistik sind jedoch im allgemeinen besonders 
die kurzen Messzeiten von Bedeutung. Hier muss man natiirlich 
das mit dem delay line clipping erreichbare Verhialtnis von Signal 
zu Untergrund mit dem des optimalen Verstiirkers bei gleicher 
Messzeit ty vergleichen. Als Resultat einer entsprechenden Rech- 
nung erhalt man in sehr guter Naherung: 


Varlte) Voptiiwi- Set, <i,2). 


Man kommt also offensichtlich bei Ionisationskammermessungen 
mit den gebrauchlichen Methoden der prinzipiell tiberhaupt még- 
lichen Messgenauigkeit schon so nahe, dass eine wesentliche Ver- 
besserung nur durch ganz andersartige Betrachtungen und Ideen 


méglich ist. Aber vielleicht hat auch diese Erkenntnis einen gewissen 
Wert. 
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